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ABSTRACT 

; A UHF clear-air Doppler radar wind profiler was purchased by the Naval PoSt- 

I graduate School and began operation in late-March of 1990. The identification of nu- 

t 

merous meteorological phenomena utilizing UHF Doppler radar wind profilers has been 
1 successfully documented; however, published work to-date has concentrated on inter- 

I pretation of hourly-consensus data. But with the onset of improved technology, it is 

I « now possible to retrieve reliable data within the hourly period. This study examines the 

quality of the six-minute data through statistical analysis and determines the overall ac- 
- curacy of the School's profiler. Specifically, the performance of the profiler is compared 

with rawinsonde data, other profiler comparison studies, as well as the contract specifi¬ 
cations. 'The study also looks at overall performance of UHF Doppler radar profiler in 
an environment characterized by a strong marine inversion. Finally the study presents 
evidence that meteorological phenomena which exhibit time scales on the order of min¬ 
utes can be observed using the six-minute profiler data. 
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I. INTRODUCTION 


Radar technology has undergone continuous refinement since it was first developed 
early in this century. The, development of the modem wind profiler is an outgrowth of 
research done in the I920's and I930's with radars designed to probe the ionosphere. 
This research was driven by the need to calculate HF radio interference, for example, 
from the velocity of the free-electrons and ions. Interest faded through the late 1930's. 
It was not until World War II that work began again with short-wave radars-hoping 
again for the chance to utilize the ionosphere successfully for communications. Al¬ 
though the initial applications of radar were for communications research and detecting 
solid objects, such as warships, it was soon discovered that radar could also detect 
hydrometeors in the atmosphere. It was not until the early 1970's that strong returns 
were first received in "clear air" (Woodman and Guillen, 1974), 

In the late 1970's the first significant work with profilers began in the United States. 
While the early wind profilers were designed to sample the stratosphere and even the 
mesosphere, today's versions concentrate on profiling the troposphere and lower 
stratosphere. The National Oceanic and Atmospheric Administration (NOAA) Envi¬ 
ronmental Research Laboratory in Boulder, Colorado began research with 50 MHz and 
400 MHz radars in the early 1980's and since then have been operating a small network 
of wnd profilers in Colorado (see inset within Fig. 1). 

For the past seven years the Denver Weather Service Forecast Office has had access 
to the profiler network and found the profiler data to be very useful in forecasting winter 
storms. Although the most experience dealing with wind profilers and their subsequent 
network capabilities has been in the Colorado area, it is expected that profilers will be 
extremely useful in other climatic regions. 

The NOAA Wave Propagation Laboratory is in the process of setting up a demon¬ 
stration network of 30 wind profilers by 1991 (locations depicted in Fig. 1 from van de 
Kamp (1988)). The first UNISYS-produced profiler was officially accepted by NOAA 
in August of 1989, located at Platteville, Colorado. As the network becomes complete, 
forecasters and researchers will have the opportunity to compute and evaluate hourly 
time series of wind profiles and derived kinematic quantities. 

The Naval Postgraduate School has purchased a 404 MHz profiler as part of the 
Marine Atmospheric Measurements Lab in Spanagel Hall. The transceiver antenna was 
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installed at Fritzchc Field, Fort Ord. The area topography and geographical position 
orihe NFS profiler are presented in Fig. 2. Figure 3 portrays details ofthe NFS profiler 
site and Hat-bed antenna array. The profiler was operational by latc-March 1990. 

The NFS profiler was tlic first UMF profiler placed on the California central coast, 
and among the first on the North American west coast. The California central coast 
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Kilometers east of Profiler Site 


Fig. 2. Topographical Map of area Surrounding NFS Profiler Site 

provides a unique environment influenced by cold-water upwelling, a persistent 
subsidence of dry air due to the eastern Pacific subtropical ridge, and the complex 
interaction of topography and the land/sea breeze circulation. Consequently, there is 
need to determine the feasibility of operating a 404 Mllz ladar profiler in an environ¬ 
ment characterized by a strong marine inversion. 

Significant work with profiler data has already been accomplished, especially by Kuo 
el at. (1987) utilizing the Colorado network and Forbes ei aL (1989 and 1990) with the 
Penn State profilers. Kuo ei al. found that retrieved temperature derived from a network 
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Fig. 3. Configuration of the NFS Profiler Site: (Source: Tycho Technology, Inc., 
Wind Profiling, 1988) 
























of VHF wind profilers was significantly more accurate than the direct temperature 
measurement from a combined satellite- and ground-based microwave radiometric sys¬ 
tem. Forbes et al. documented how wind profiler data together with conventional 
National Weather Service rawinsonde data proved quite helpful in understanding the 
cause of subtle weather systems seen in satellite imagery. Forbes et al. also examined 
the utility of profiler-derived temperature gradients, temperature advections and lapse 
rate gradient information. Both the Colorado and Pennsylvania network of profilers 
have been operating for nearly a decade and, subsequently, numerous works, beyond 
which can be cited in this study, have analyzed important aspects of their respective local 
meteorology, including frontal passages and convective storms. 

Published studies to-date have concentrated on interpretation of hourly-consensus 
data. Forbes et al. (1990) discusses, although only briefly, interhour fluctuations de¬ 
tected with the Penn State VHF profiler, but their emphasis is on analyzing significant 
mesoscale weather. The primary objective of this study, in contrast, is examination of 
the potential for a clear-air wnd profiling Doppler radar to observe small-scale, high- 
frequency meteorological cvcnls-in other words, information within the time scale of 
one hour, down to time scales of six nunutes in the case of the NPS radar. To meet this 
objective it will be necessarj’ to first ensure data from the individual profiles are as reli¬ 
able as possible, i.e., throw out the "bad" data. This data quality process will hopefully 
provide an improved capability to obtain good representative data beyond that which 
the manufacturer's software provides. It is also hoped that by improving the individual 
profiles that the hourly-consensus will be improved. 

A secondary objective of this study is to document the performance, including ac¬ 
curacy, of the NPS profiler operating in an environment characterized by a strong ma¬ 
rine inversion. Determination of system accuracy is important in this study in order to 
lend insight into the thermodynamic relationships presented later. The study will also 
compare the data accuracy to other profiler and rawinsonde studies. 

Finally, this study hopefully will expose many students and instructors to the capa¬ 
bility of a UHF radar profiler, in addition to providing future students the necessary 
groundwork to expand and analyze specific meteorological phenomena with high- 
temporal resolution and on a real-time basis. 

This study documents research performed during the period of initial operation, 30 
March to 05 May 1990. The period was selected for several reasons. First, the profiler 
was scheduled to be part of Tropical Cyclone Motion experiment (TCM-90) in southeast 
Asia, which meant disassembly by no later than June 1990. Second, a transmitter am- 
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plifier failed one month prior to scheduled disassembly, and the amplifier subsequently 
had to be removed and sent back to the manufacturer. The period left (a little more 
than a month) was further reduced in length by spurious software failures and power 
outages. 

Several other factors limited the ability to conduct a detailed meteorological inves¬ 
tigation. The initial operating period was never meant to be analyzed in such detail as 
this study is directed. Therefore, the ancillary data-gathering systems such as SODAR 
and high-frequency pressure and humidity recording devices, were never arranged. 
Rawinsondes were the only equipment configured for this experiment, but at the last 
moment the primary rawinsonde system allotted for comparison studies failed. A sec¬ 
ondary unit was available but could be used only intermittently, as it had to be scheduled 
around cruise course requirements. In the end, only nine successful rawinsonde 
launchings were available, and five of those were from onboard the student cruise boat 
■'•st outside Monterey Bay. 

The equipment configuration used for this study is described in Appendix A. The 
Radio-Acoustic Sounding System (RASS) modification for the wind profiler was not 
installed for this period, and the software for Q computations was not yet available 
from the manufacturer. 
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II. BACKGROUND 


A. THEORY OF PROFILER OPERATION 

The wind profiler, like all radars, detects the transmission or reflection of 
electromagnetic radiation at certain frequencies, or wavelengths. A phenomenon known 
as Bragg scatter is critical to receiving the necessary backscattered signal strength in 
clear air. Bragg scatter comes about in two ways: (I) through the scatter from a surface 
or layer with a wave speed which is resonant to the cosine of the incident angle of the 
transmitted radar pulse, or (2) through the presence of multiple refractive layers spaced 
at the scales of turbulent energy as is in the case of UHF profilerS“404 MHz falls in the 
Ultra High Frequency (UHF) range. Although all scales of turbulence scatter energy, 
it is not until the turbulence scales approach 1/2 the transmitted wavelength, or 37 cm 
in the case of the 404 MHz profilers, that the scattered energy waves interact construc¬ 
tively. 

Profilers depend upon the scattering of electromagnetic energy by minor irregulari¬ 
ties in the index of refraction of air. Energy is first radiated at 404 MHz from a coaxial 
colincar antenna (Fig. 4), essentially several cylindrical patterns radiating from inter¬ 
leaved coaxial cabling. Above approximately 200 m, beyond the near-field, these indi¬ 
vidual cylindrical patterns begin to interact constructively to form the radio wave front. 
The antenna is constructed so most of the energy is embodied within a five degree cir¬ 
cular beam, or cone; although several side lobt arc present as illustrated by the insert 
within Fig. 4. This wave front is not impeded intil turbulence on the correct scale is 
encountered, as shown schematically in Fig. 4. Turbulence produces strong small-scale 
variability in three dimensions, carrying moisture and temperature discontinuities verti¬ 
cally (and horizontally) adjacent to other air parcels. Atmospheric turbulence is gener¬ 
ated at scales of tens to hundreds of meters by wnd shear and convection. These large 
scale eddies break up into smaller eddies, which break up into still smaller ones, and so 
on until the turbulent energy is dissipated at very small scales by viscous heating. This 
cascading of kinetic energy occurs partially within a wavenumber space (ranging from 
only a few centimeters to hundreds of meters) referred to as the "inertial sub-range". It 
is somewhere within this inertial sub-range where turbulent scales of the multiple of 
one-half the transmitted wavelength are encountered. When an electromagnetic wave 
front encounters these areas of variability, the wave front is refracted slightly and a small 
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amount of energy is scattered in all directions, as represented by the dotted wave fronts 
emanating from point 'AA' within turbulent layer 'A' in Fig. 4.1 Very high above the 
troposphere the density of air is such that it cannot support small scale turbulence, 
hence the first theoretical limitation on the profiler is a maximum altitude that can sup¬ 
port turbulence on the scale of 37 cm-a height which varies but is estimated to be be¬ 
tween 18-20 km for the West coast. 

The degree of turbulence and the variations in refractivity index over short distances 
affect the energy backscattered (the power returned). The refractive index for the 
troposphere at UHF frequencies (Balsley and Cage, 1980) is 

/i=I-f77.6x 10*^ Y+3.73x10“'-^, (1) 

where e is the partial pressure of water vapor in mb (humidity), p is the atmospheric 
pressure in mb, and T is the absolute temperature in K. A refractive index structure 
function Z),(r) described by Gossard and Strauch (1983) is used to quantify these vari¬ 
ations over short distances as 

( 2 ) 

where x and j? -h v are the vectors originating at the profiler antenna and terminating at 
the centers of two air parcels (depicted as circles) separated by a distance and direction 
f, represented in Fig. 5. Within the inertial sub-range of homogeneous, isotropic tur¬ 
bulence, Gossard and Strauch (1983) state that 

Zl,(7) oc (3) 

To make an equality, a standard parameter Q , or the "refractive index structure pa¬ 
rameter", was conceived. Q is a constant of proportionality describing the refractive 
index differences between two parcels of air situated adjacent to one another; so now 
Eq. (3) becomes 

D,(r)=Cy<^. (4) 


1 There is evidence suggested by Doviak and Zmic (1984) and Waterman et al. (1985) that at 
certain wavelengths, the backscattered energy is not isotropic, but rather anistropic. There is on¬ 
going research to resolve this phenomenon. 
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Fig. 4. Illustration of the Generation, Beam Pattern, and Scattering of Profiler 
Emiltcd Signals 

Q Is dependent on the temperature and humidity, as well as the pressure, as given by 
the relationship from Peterson (1988) 


C^= 79 X 10~*[-^][Cf + .113C7-2+ 3.2 ^ 







Fig. 5. Illustration of Refractive Index Structure Sampled by the Profiler for a 
Singic-Scattcrcr 

where Q- is the “temperature structure function”, Cq is the “humidity structure function”, 
and Cj(, is a combined structure function. In the inertial sub-range, the radar 
"rcllcctivity” (tj) is proportional to Q, and is given by Battan (1973) 

1/ == 0.38CJr/l“‘'^ (6) 

where A is the radar operating wavelength. 

Radar reflectivity can now be expressed as returned power (P,), derived as Uq. (B-22) 
in Appendix B 


= (7) 

R 

where K is a "radar constant" (i.e., radar efficiency and mathematical constant), P, is the 
power transmitted, t is the pulse duration, A, is the clTcctive antenna area, R is the dis¬ 
tance from the antenna to the turbulent eddy, and tj is the radar reflectivity from Cq. (6). 
R is determined by the formula R = ct/2, i.e., the lime it takes for backscattered energy 
from the respective sampling heights to return to the antenna; c is the propagation speed 
of energy at radio frequencies and t is the time of transmission subtracted from the time 
of receipt. 
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B. PROFILER DATA PROCESSING 

The wind speed and direction can be derived trigonometrically from the radial ve¬ 
locity of the profiler's individual transmitted beams, of which there are three-one verti¬ 
cal and two off-vertical. Figure 4 graphically depicts how these beams are formed; the 
off-vertical beams are phased by delaying the transmitted RF pulse to certain portions 
of the array. 'The radial velocity is derived from the Doppler frequency shift created by 
turbulent eddies traveling along with the mean wind. In the presence of turbulence or 
other velocity inhomogeneities (refractive index fluctuations) within the beam, the return 
signal spectrum may vary appreciably from one pulse to the next. In order to measure 
the prevailing flow velocity (by averaging out microscale fluctuations), many samples 
(about 1400) are composited over a period of integration (i.c., sample period). Even if 
there were no velocity variations within the beam, the returned signal would be at an 
assortment of frequencies, since the transmitted signal is not of uniform frequency. (In 
addition to the 404 MHz frequency, the transmitter also leaks low-intensity "tails" at 
slightly higher and lower frequencies. As the pulse length decreases the bandwidth in¬ 
creases.) 

The period of integration is 58 seconds in the low mode and about 55 seconds in the 
high mode. The result at the end of the integration period is the most fundaniental 
profiler quantity, a vertical profile of 36 wind velocity spectra. A representative example 
of the velocity-power spectra is given in Fig, 6. There are six modes of data: two off- 
vertical and one vertical mode for each the high power and low power configurations. 
In order to sample the entire troposphere and lower stratosphere, low and high mode 
profiles are measured sequentially, requiring a total time of approximately six minutes. 
Thus, 10 sets of tropospheric winds can oc obtained p3r hour. 

To get data, a radio wavefront or RF pulse must first be transmitted. After the en¬ 
ergy is radiated by the antenna system, an electronic circulator switch is thrown to re¬ 
connect the antenna to the receiving equipment. There are 100-150 microseconds 
between each pulse and it is during this interval that the receiver samples the returned 
signal at different "gates" or vertical range bands as illustrated at the top of Fig. 7. 
Then, for each gate, several thousand samples are gathered and time domain averaged 
by the processor, as shown in the second and third functional steps in Fig. 7. 

Time domain averaging (TDA) consists of taking the mean of several consecutive sam¬ 
ples in time. Figure 8 functionally illustrates how the received energy is broken into 
quadrature and in-phase components, and the subsequent time-averaging of each com¬ 
ponent. This time domain averaging is applied to the power spectra of the 36 available 
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Fig. 6. Typical Vertical Profile of 36 Range Gate Frequency-Power Spectra: (As 
displayed on the Keyboard Display Unit (KDU).) In the upper-left hand 
corner of each gate is the center height of the sampled volume in meters. 
Beneath the height parameter is the resolved radial velocity in m s'\ also 
indicated by the vertical bar centered within the dominant peak. In the 
upper-right of each spectrum is listed the .signal-to-noise and below this 
parameter is the ma.ximum pow.r encountered within the spectrum, both 
in dli. 


gates, hence a TDA matrix can be formed as illustrated in Fig. 8, which consists of 256 
samples in each of the 36 gatcs"256 samples are needed for the fast fourlcr transform 
(FFT) calculations later. The entire time the processor takes to collect a full matrix for 
all range gates is approximately three seconds in the low mode. 

'fo this point in the data processing, a typical power spectrum in relation to time is 
available for each gate. But what is desired is the dominant frequency, therefore fre¬ 
quency domain averaging or spectral averaging subsequently ensues. It consists of cal¬ 
culating a mean spectrum from several consecutive spectra after the FFf computations 
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Schematic Diagram of the Primary Signal Processing Steps for the Profiler 
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Fig. 8. Funtional Diagram of the Time Domain Averaging (TDA) Process 


have been performed. Now the time required to complete a spectrum in each of the six 
modes can be determined: 


Time_per_spectrum = PRPx TDA x 256, 


where PRP is the pulse repetition period of the radar and TDA is the number of time 
domain averaging matrices. Finally, the number of spectra utilized in the averaging 
process of each mode is adjusted, so the total time per spectrum is approximately 60 
seconds, and, consequently, the total lime per profile is then six minutes, which matches 
the National Oceanic and Atmospheric Administration reporting requirements for their 
future profiling network. The number of spectra utilized for a single profile during the 
first month of operation is listed in Table 1. 

At this point, a typical power spectra in the frequency domain has been processed 
for each mode at each gate from which wind velocities must be determined. Briefiy, a 
wind profiler is a Doppler radar that measures the shift in frequency created by moving 
targets. In the case of the wind profiler, the targets are turbulent eddies moving with the 
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Table 1. PROFILER PARAMETERS 


Mode 

High 

Low 

No. of 
TDA Ma¬ 
trices 

40 

120 

PRP 

162ps 

100/rs 

Tinte per 
Spectrum 

1.65s 

3.072s 

No. of 
Spectra 

34 

19 


mean wind. Tims by determining the Doppler sliift of the spectral moments, the wind 
toward and away from the radar can be determined. To accomplish this, a 
manufacturer's-supplicd algorithm computes the Doppler spectral moments of the 
spectral peak a.ssociated with a local maximum in the frequency spectrum. Inputs to this 
algorithm are the velocity (or frequency) versus power spectrum (with noise and ground 
clutter removed) and the index of the local maximum determined to be signal. Outputs 
of this algorithm arc the three spectral moments: 

1) Signal Power (estimated total power of the received signal), 

2) Mean radial velocity (first moment of the spectral peak), and 

3) Radial velocity variance (second initial moment, or width, of the spectral 
peak). 

I his routine is specified for each mode and for each gate. The Doppler shift of the peak 
(f„) can be converted into a radial velocity component of the wind, by rc-afranging the 
following basic equation from Battan (1979): 



'fhe central frequency (as shown in Fig. 7) corresponds to the transmitted frequency. 
If a "target" is moving toward the radar, the frequencies will be increqsed; conversely, if 
it is moving aw'ay, the frequency will decrease. The Doppler return and the width of the 
peak at each range gate arc measures of the degree of turbulence in the sample volume. 
'I he second peak within the sample spectra shown in Fig. 7 is an indication of detected 
precipitation. (Precipitation would be verified by on-site observations.) 
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Radial velccity components are next decomposed into easterly and northerly com¬ 
ponents of the Cartesian coordinate system. This is accomplished by the vertical and 
two off-vertical beams. The vertical beam pattern measures the actual w-component of 
the wind. With the w-component known, the u- and v-components can then be derived 
by trigonometry, as shown in Appendix C. Within Appendix C, and 6j (the angles 
of the offset beams) are most commonly set at 15* from vertical and tf) (the angle in the 
horizontal plane between the two off-vertical beams) is set at 90* for case of trigono¬ 
metric calculations. 

Finally, the three moments are consensus-averaged, based upon data collected for 
a one-hour averaging period. The consensus technique is detailed in Appendix D. As 
suggested by the name, not all the data are averaged, rather only those that mutually 
agree within a threshold. A 2 m jt* window, or threshold, is used by the manufacturer 
for the radial velocities. If the consensus group does not have at least four members, 
then the average value is reported as missing. This subset is then used to calculate the 
average values of all three moments. 

C. PROCESSING OF RAWINSONDE DATA 

The rawnsonde signal is sampled by the ground system every ten seconds. The wind 
data from 24 consecutive ten-second samples arc then averaged-this provides smoothing 
of the Initial data to remove noise (see Fig. 9). Consequently, all sharp features are ef¬ 
fectively removed. A running average is subsequently maintained with data output every 
40 m, based on a four mr' balloon ascent rate. A detailed presentation of the scheme 
for deriving rawinsondc data carmit be given within this study because the material is 
proprietary with Vaisala Corporation. (Winds in the lowest 200 m are handled in a quite 
different fashion, due to the obvious lack of data required for a running average; this 
different technique is proprietary information but is not important for this study since 
there are no profiler wind values at these levels.) 

The rawinsonde data then represent approximately a 1000 m (actually 960m) thick 
volume (with 40 m gate spacing between subsequent data), whereas the wind profiler 
has 250 m resolution, or interval spacing, in the low mode (with 250 m gate spacing) and 
1000 m resolution in the high mode (with 500 m gate spacing), hence different values as 
illustrated in Fig. 9. Therefore, a typical height of rawinsonde measurements centered 
within the sampling profiler gate was required for comparison. U and V-velocity 
rawinsonde components were subsequently averaged through a column centered around 
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the desired profiler gate height. Likewise, rawinsonde temperature and humidity were 
linearly averaged with respect to height around the profiler gate heights. 

D. KNOWN SOURCES OF PROFILER ERROR 

There are many potential sources of noise in the profiler radial velocities. Of most 
significance for this comparison study, the profiler and rawinsonde sample at different 
locations. The wind profiler provides simultaneous measurements over all heights at a 
fixed location. It samples large volumes of atmosphere (up to 1 km deep) and assigns 
an integrated value to the center height of the sampled column. The rawinsonde, on the 
other hand, records measurements at different, discreet heights and at different times (a 
typical ascent will take up to 1-1/2 hours and may cover many kilometers lateraUy). As 
a result, the rawinsonde and the wind profiler often obtain measurements at widely sep¬ 
arated locations where the winds may be different. Also, the winds over the fixed profiler 
location are not uniform during the rawinsonde's 1-1/2 hour ascent; therefore, variability 
of the atmosphere can make the interpretation of any comparisons difficult at times. 
Furthermore, since these instruments sample in a functionally different manner, com¬ 
parisons of their measurements cannot provide unambiguous information on the accu¬ 
racy of cither instrument. 

Another large source of error is electromagnetic noise present within the pulse vol¬ 
ume as a result of nearby power lines or machinery, as well as other contributions fror.i 
side lobes (power “leaking' from the antenna at angles other than the desired angle >. 
The Federal Aviation Administration began operating a traffic control radar within 60 
m of the Ft Ord wind profiler site in April of 1990. A'ihough the operating frequency 
docs not appear to affect profiler operations, the movement of the antenna provides a 
strong enough Doppler shift of the profiler transmitted frequency to provide a significant 
return in the side lobes of the receiver entenna and contaminate the lowest gate. There 
appear to be many forms of this contamination such as the inclusion of a sinusoidal 
wave pattern in the lowest four gates as represented in Fig. 10. A more common con¬ 
tamination pattern is the occurrence of narrow peaks at approximately 1ms-' (not 
shown). 

Noisy wind vectors can result from an imbalance of the individual beam sensitivities, 
since one beam may .detect the weak atmospheric turbulent signal while the other may 
intermittently resolve the background noise as the dominant signal source. In the event 
that the atmosphere is quite dry and the associated profiler backscatter w'eak, there is a 
fairly strong persistent background noise frequently present and the “insensitive' beam 
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Fig. 9. Functional Ran insonde and Profiler Gate Spacing 


may 'lock in on" the background frequency. The result is a high consensus nuntber but 
an unrealistic wind. 

Closely related to this latter error is the susceptibility of each beam to be less sensi¬ 
tive to cross-winds, or those prevailing winds which cross perpendicular to the direction 
which the beam is formed. Hence, the accuracy of determining these cross-wind veloci¬ 
ties is reduced. There is consistent evidence of this bias throughout the course of this 
study, and analysis of a larger data set would probably define the bias discreetly. 

When the radar switches from transmit to receive mode, there is electronic noise in 
the system that takes a few microseconds to dissipate. Any signals returned by the at¬ 
mosphere during this brief period are masked by noise. Because the first signals to re¬ 
turn come from low in the atmosphere, the lowest sampling gate (446 m) is most 
susceptible to this problem. No evidence of this elTcct was discovered through the 
course of this study. 

The profiler operating frequency and sampling interval determine the Nyquist ve¬ 
locity, the maximum radial velocity that can be measured. The well-known Nyquist 
sampling theorem states that the highest frequency to be resolved in the signal must be 
sampled at least two times per period. If the signal contains energy at frequencies higher 
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Fig. 10. Sample of FA A Traffic Control Radar Contamination of NFS 
Profiler: Evident in the four lowest gates. (See Fig. 6 for format de¬ 
scription.) 


than the Nyquist frequency, the Doppler signal can be aliased by higher frequencies. 
'I’his means that the energy is "folded back" to frequencies less than the Nyquist fre¬ 
quency (McGlllcn and Cooper, 1974). For the NFS profiler the Nyquist velocity is 
about 15 ms~' (30 kt) in the low mode and 60 ms*‘ (120 kts) in the high mode. Radial 
velocities greater than the Nyquist velocity arc "folded" as illustrated in Fig. 11 at 11697 
m. Although velocity folding occurred in two spectra within Fig. 11, fortunately it was 
not folded into the spectral region "sensitive" to the algorithm. There were no known 
occurrences of this type error during the study period. 

Contributions from persistent vertical velocity, such as in the vicinity of sloping 
terrain or the closed cells of the land/sea breeze circulation, may cause appreciable con¬ 
tamination of the "horizontal" velocity measured by the off-zenith beams, especially 
when the horizontal wind speed is low. Since the "horizontal" beams are 15® from 
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Fig. II. Sample of Velocity Folding in the NFS ProFilen Seen at the 11697 m 
level. (See Fig. 6 for format description.) 


zenith, 75 percent of the oflT-zenith radial velocity is sensed as vertical velocity. Thus the 
horizontal velocity required to produce a 1 mj~' radial component is about 4 nts’K 
Synoptic-scale vertical velocities (~ O.I ms*') appear not to pose a problem, but 
smallcr-scalc systems which may contain vertical velocities as great as several meters per 
second arc clearly a potential source of horizontal beam contamination. 

Turbulence accompanied by significant backscatter above a height which the profiler 
can "sec" introduces aliasing into the data. The ma.ximum unambiguous radar height is 
given by Battan (1974) as 


R, 


ciPRP) 


mix 


where c is the propagation speed of radio energy. PRP as set in the low mode of the 
NPS profiler is 100 fis, which sets in the low mode at 15,000 m. (/?„„ in the high 
mode is 24,300 m.) Therefore, if the transmitted signal is returned from a source higher 
than 15,000 m, the radar processor while in the low mode will sec this signal as a possible 
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early return in the lower gates of the next sampled pulse. Consequently, ionospheric and 
cosmic noise can also b,e problems. In this case the radar will receive noise in any of the 
gates, when realistically the signal evolves quite a bit higher beyond that which the radar 
can measure. 

Precipitation and other hydrometeors can also scatter energy back towards the an¬ 
tenna. Typical meteorological radars operate at wavelengths on the order of 5 to 10 cm, 
whereas wind profilers operate typically near wavelengths of 74 cm. Hydrometeors are 
some 10* times less reflective to wind profilers compared to conventional weather radars. 
This is indicated by the reflectivity equation given by Gossard and Strauch (1983) for 
single particle scattering: 


{const)t^ 



where d is the particle diameter. Precipitation, when detected, can be a problem in some 
circumstances. Clearly, if there is interest in the vertical motion of the air, precipitation 
presents a problem, as the vertical velocity measured is that of the precipitation, and not 
the air. In some cases, spectra from the vertical beam will show a double peak: one 
corresponding to the clear air and the other precipitation. This is illustrated in Fig. 12. 

Rawinsondc operation colocated with the profiler creates a false wind velocity value. 
Figure 13 presents two spectra from the Keyboard Display Unit. Balloon interference 
is highlighted at approximately 6000 m in the top frame and 11000 m in the bottom 
frame. The lime diflerence between these two frames is 25 minutes. Again, assuming a 
A rnsr' balloon ascent rate, this is precisely the height the rawinsondc should have at¬ 
tained. The detection of the balloon appears to be due to the difTereni index of refrac¬ 
tion value characteristic of the gas within the balloon vice the actual reflection from the 
non-metallic surface of the balloon. The radio signal of 402 MHz utilized by the 
rawinsondc for passing data to the ground system does not appear to interfere with 
profiler operation. 

It must also be recognized that the vertical beam and the two off-zenith beams do 
not sample the same volume of air, in either space or time. There may be up to a 3.2 
km difference at an altitude of 12 km. This is especially troublesome in the presence of 
small-scale mesoscale disturbances. 
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Fig. 12. Example of a Double Peak in NFS Profiler Spectra causediiy Rain; ev¬ 
ident clearly at the 914 and 1214 m levels, as well as other disparities ex¬ 
tending up to 4000 m. (Sec Fig, 6 for format description.) 
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in. DATA ANALYSIS 


A. COMPARISON OF MANUFACTURER'S CONSENSUS TO RAWINSONDE 
DATA 

Prior to analyzing dynamic, kinematic and thermodynamic relationships derived 
from the NPS profiler winds, it is necessaiy to first determine the accuracy of the profiler 
wind data. This information provides the NPS Marine Atmospheric Measurements Lab 
a benchmark of how the NPS profiler performs in relation to other profilers as well as 
verifying requirements within the profiler procurement contract. Both the wind profiler 
and rawinsonde system provide vertical profiles of wind measurements; therefore, like 
many previous studies, the accuracy of the profiler will be gaged by the rawinsonde data 
taken as truth data. A statistical comparison between profiler and rawinsonde wind data 
was made and used two approaches: (1) each horizontal wind component was correlated 
between the profiler and rawinsonde as determined by a correlation coefficient, and (2) 
mean and standard deviations for the difference in each horizontal wind component as 
measured by the two instruments was calculated. Because such a small rawinsonde data 
set was available during the first month of operation, a more detailed look into why 
potentially "bad" profiler data was not in agreement with the rawinsonde data was done 
manually, providing, a third, more subjective approach. 

Nine individual profiles are compared to assess the quality of the manufacturer's 
consensus winds measured by the profiler. Figure 14 is a representative comparison of 
the vertical profiles of wind speed and direction measurements for 20 April 1990 made 
by the two instruments. The associated complete rawinsonde sounding is presented in 
Fig. 15. Additional comparisons for the remaining eight soundings (Figs. E-I through 
E-16) as w'ell as the specific positions of the five at-sea rawinsonde launchings (Fig. 
E-17) arc presented in Appendix E. In each case the rawinsonde profile is flanked by its 
associated Tycho hourly-consensus profiles, approximately 2 hours on either side of the 
rawinsonde launch lime. Subjectively, the profiles correlate very well, at least where 
consensus in the radar data is achieved. The wind direction correlates very well. Below 
5 km the wind speeds are generally within 5 kts of each other, while above 5 km the wind 
speeds are accurate to within 10-15 kts. The departure between individual profiler and 
rawinsonde wind barbs in the lowest 2000 m seen from the remaining launch times in 
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Appendix E could potentially be due to the discrepancy in instrument locations as well 
as the turbulence detected by the radar profiler in the presence of the marine inversion. 

A comparison (rawinsonde vs profiler) of the individual horizontal wind components 
versus height are illustrated in Fig. 16 and 17 for a representative low mode profile from 
0000 UTC 20 April 1990 and Fig. 18 and 19 for a typical high mode profile from 0700 
UTC 03 May 1990. There is generally very good agreement below 5-6 km in the low 
mode and below 8-9 km in the high mode. There is a notable bias in the v-velocity 
component in the high mode, which appears to be due to the orientation of the prevail¬ 
ing winds in relation to the profiler antenna, which in this case were from the north. 

In the low mode the maximum height obtained appears to correlate with the fall-off 
in humidity (evident in Fig. 15 where the atmosphere becomes dry above 4000 m); but 
is never greater than approximately 8 km even when humid layers are present aloft. This 
is due to the weakened signal from geometric spreading. (A humid layer aloft is defined 
in this paper as a layer with greater than 40% relative humidity above 4000 m.) This 
limitation is expected since a .74 cm radio signal is four times more sensitive to humidity 
fluctuations associated with turbulence than to temperature fluctuations, at least in the 
lower atmosphere (Cage 1990). This is illustrated in Fig. 20, by a deep humid layer 
which exists to 5 km. This is reflected in both the U and V-velocity component winds 
in Fig. 18 and Fig. 19. 

As additional proof of the dominance of humidity in the lower atmosphere, a nearby 
sounding was compared with the raw spectra. Fig, 20 is the 1200 UTC 06 April 1990 
Oakland sounding, 10 hours prior to the low-power mode profiler spectra in Fig. 21 
(0000 UTC 07 April 1990 was not available). Typical conditions for this period of study, 
as seen in Fig. 21, are first a layer of very strong returns between 500-1400 m, the ma¬ 
rine thermal inversion; which is followed by a rapid fall-off of signal return, until ap¬ 
proximately 3000 m is reached where the signal levels begin to increase due to humidity. 
This corresponds very well with the moisture layers in the sounding. This is then fol¬ 
lowed by a fall-off between 4000-5000 m and a local maximum of -2 dB S/N reached at 
5000-5500 m. In summary, a very good subjective correlation between the rawinsonde 
sounding and profiler spectra is evident. 

The maximum height attained in the high mode is about 12 km (significantly less 
than the theoretical 24 km which the electronics allow), again due to weakened signal 
return from geometrical spreading. But, what is additionally seen in the U and V com¬ 
ponents is a 2-3 km gap in consensus data just below the tropopause (see Fig. 18 and 
Fig. 19). The tropopause is identified by a constant temperature profile seen in the 


25 




rawinsonde data. It appears at this time, that this gap is centered around the maximum 
of the "jet stream", or jet stream core if overhead. This is the most consistent feature 
seen in the data above 8 km. Apparently, the wind shear throughout this region is 
minimal, thus, reduced turbulence in a relatively dry region provides little signal return. 
However, above the jet maximum, winds begin to decrease and radar returns are de¬ 
tectable above noise. This too is seen in the raw spectra in Fig. 22, where signal 
strength falli off and then increases above -10 dB at the 12197 m level. 

To quantify the data accuracy, a statistical comparison was made by computing the 
correlation cocITicient, calculating the standard deviation, and graphically portraying the 
data "pairs" on scatter plots. The comparisons were based on unedited consensus data 
from the manufacturer's algorithm. Figures 23 and 24 are scatter plots of component 
rawinsonde measurements versus wind profiler measurements for the nine combined 
launchings. Each marker corresponds to a measurement made at about the same hour 
and height for both instruments. If the two measurements are equal, then the marker 
will lie on the diagonal. If the marker is positioned off this diagonal the two measure¬ 
ments arc different. Consistent markers to one side or the other indicate a bias. 

In the low mode, the density of data around the scatter plot diagonal appears good 
below 7-8 msr' (15-17kts). Winds above this value are generally associated with higher 
heights; and as Fig. 23b illustrates, the "bad" values (the values that do not fall closely 
along the diagonal) generally lie above 6 km as represented by the markers O and A. 
Similarly, in the high mode data depicted in Fig. 24b there appears to be an appreciable 
density of data along the diagonal, uith any "bad" data related to speeds of 8 m j-* 
(17kts) above 8 km. The depiction of "good" data at higher speeds indicate that height 
is the primary contributor to the presence of "bad" data. There is a bias tendency in 
v-velocity components which is more noticeable in the low mode data. This cannot be 
explained but may be due to such a small data set, or due to the orientation of the pre¬ 
vailing winds relative to the profiler fixed beams. 

Table 2 presents the statistical results broken out by atmospheric layers. The values 
listed under the column "Mean", represent the mean error away from the mean wind 
component, while "Standard Deviation", is that value taken around the mean wind 
component, within one standard deviation or 67% of the data. The "Correlation Coef¬ 
ficient" values were calculated by the Pearson correlation coefficient relating profiler to 
rawinsonde data. 'The guidelines this study has set for "good" data is a correlation co¬ 
efficient above .93 and a standard deviation less than 2.5 m s~' , which are taken from 
the findings by Weber and Wuertz (1990). However, the performance of the radar in 
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relation to the rawinsonde system never quite met these criteria. The table reveals good 
correlation in the low mode below 6 km in both components. Good correlation in the 
high mode is seen below 8 km. There is poor correlation in the lowest 2-3 km under the 
presence of the inversion layer, i.e., there is turbulence created by the wind shear, the 
inversion jet, and/or convective thermal cells. The statistical results appear to reinforce 
the conclusions made earlier by the subjective analysis. 

An accuracy determination would not be complete without the inclusion of the site 
survey, which establishes exactly where the antenna is oriented. For consideration, the 
elevation accuracy of the antenna base height is ± 1 m, while the azimuthal accuracy is 
± 3' of arc. The results and the survey site layout are presented in Appendix F. The 
accuracy to which the profiler array was established is considered negligible for this 
study, but is presented for completeness and for the record to accommodate future 
studies. 

The NFS wind profiler appears from the statistical analysis to be accurate to within 
only 3-4 m r-' in the lower atmosphere (or more specifically within the moist layer below 
6 km). This error is slightly greater than the 2.5 m r-' standard deviation computed by 
Weber and Wuertz (1990). The profiler error is significantly greater than the 1 mr' 
accuracy specification listed in the procurement contract. Meanwhile, Strauch et al 
(1987) found a standard deviation of about 1.3 mj-' in clear air and Wuertz ei al. (1988) 
found a standard deviation of about 2-4 m r' in precipitation. Thomson and Williams 
(1990) reported amazing accuracies in profiler data of ±l msrK The studies by Weber 
and Wuertz, as well as the other studies presented, analyzed several thousand pairs of 
rawinsonde and profiler data, while this study only was able to look at a few hundred. 
A larger data set would be expected to reveal a lower standard deviation based on the 
small mean error values listed in Table 2, for heights below 6 km. 
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Sonde High Mode Comparison Data: 90 04 19 at 23 23 GMT 
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Fig. 14. Pronier*Ra>unsonde Comparison Winds for 0100 UTC 20 April 1990 


28 


Profiler Profiler Profiler Profiler 





Lou nch time : 90 04 19 al 23 23 GMT 


Altitude (meters) 




o 

o 

O 

o 

o 

o 

o 

o 

o 

o 

o 

o 

O 

o 

o 

o 

o 

o 

o 


o 

o 

O 

o 

o 

o 

o 

o 

o 


03 




o 

CO 



IN 



(sjaiaiu) apnijilV 


15. Rawinsonde Sounding for 0100 UTC 20 April 1990 
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Fig. 16. Typical Comparison Plot of Ran insonde vs Profiler U-Velocity Component 
(Profiler in Low-Mode): Rawinsondc winds are indicated by the o 
marker, while the profiler-derived winds are indicated by the Hh marker. 
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Fig. 17. Typical Comparison Plot of Ran insonde vs Profiler V-Velocity Component 
(Profiler in Lon-Mode): Rawinsondc winds are indicated by the o 
marker, while the profiler-derived winds are indicated by the + marker. 
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Fig. 18. Typical Comparison Plot of Ra>viiisonde vs Profiler U-Velocity Component 
(Profiler in High-Mode): Rawinsondc winds are indicated by the o 
marker, while the profiler-derived winds arc indicated by the + marker. 
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Fig. 19. Typical Comparison Plot of Rawiiisonde vs Profiler V-Velocity Component 
(Profiler In High-Mode): Rawinsondc winds arc indicated by the o 
marker, while the profiler-derived winds are indicated by the + marker. 
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Fig. 20. Oakland (72493) Sounding for 1200 UTC 06 April 1990 
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Fig. 21. Example of Lo>v Mode Spectra: There is increased signal strength in the 
5196 and 5446 m levels due to humidity. 
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Fig. 22. Sample Lotv Mode Spectra Depicting Jet Stream Structure: Evident is 
an increase of signal strength at the 12197 m level. 
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ig. 23. Scatter Plot of Lo>v-Mode Profiler vs Rawinsonde Velocity: Data taken 
from all nine rawinsonde launches. The marker x represents those winds 
derived below 2 km, o from 2 to 4 km, + from 4 to 6 km, A from 6 to 8 
km, and O any winds above 8 km. 



Fig, 24. Scatter Plot of High-Mode Profiler vs RaMinsonde Velocity: Data taken 
from all nine rawinsonde launches. The marker x represents those winds 
derived below 4 km, o from 4 to 6 km, + from 6 to 8 km, A from 8 to 10 
km, and O any winds above 10 km. 
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Table 2. STATISTICAL ANALYSIS OF DIFFERENCES BETWEEN 
RAWINSONDE AND PROFILER-DERIVED VELOCITIES FOR ALL 
NINE RAWINSONDE LAUNCHES 



wm 

Layer 

mmSmM 

Standard 

Deviation 

Correlation 

Coefficient 



All Heights 
Combined 

2.49 

10.31 

.490 



<2 km 

0.083 

3.4 

.791 

Low 

U Velocity 

2-4 km 

-.828 

4.66 

.815 

4-6 km 

-.899 

8.4 

.783 



<6 km 

-.553 

5.78 

.793 



6-8 km 

-4.59 

14.09 

.285 



> 8 km 

-13.0 

18.68 

.077 



All Heights 
Combined 

7.54 

14.31 

.415 



<2 km 

2.06 

4.04 

.634 

Low 

V Velocity 

2-4 km 

1.67 

4.31 

.757 

4-6 km 

3.10 

8.01 

.714 



<6 km 

2.22 

5.65 

.740 



6-8 km 

18.98 

19.25 

0.01 



>8 km 

26.0 

20.34 

-.048 



All Heights 
Combined 

-4.94 

12.05 

.439 



<4 km 

-.996 

3.275 . 

.912 

High 

U Velocity 

4-6 km 

-1.04 

3.091 

.970 

6-8 km 

0.627 

13.78 

.516 


‘ 

<8 km 

-.717 

6.44 

.805 



8-10 km 

-15.36 

21.43 

.063 



> 10 km 

-16.35 

14.69 

.647 



All Heights 
Combined 

4.88 

12.83 

.361 



<4 km 

1.985 

3.817 

.742 

High 

V Velocity 

4-6 km 

0.547 

3.916 

.914 

6-8 km 

4.29 

10.77 

.611 



<8 km 

1.99 

5.79 

.783 



8-10 km 

16.88 

29.61 

-.609 



> 10 km 

11.59 

18.96 

.217 
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B. INVESTIGATION OF VARIABLES AFFECTING DATA QUALITY 

Given that there are numerous sources of error in the profiler operatioii that influ¬ 
ence the six-minute profiles and hourly-consensus, it may be possible to eliminate bad 
six-minute data in order to subsequently improve overall performance. To assess the 
six-minute data, several parameters related to the basic moments calculated for each gate 
are examined: Signal-to-Noise power Ratio (S/N), beam power differential, and velocity 
variance (or width). Although one of the most effective data quality control techniques 
is a running comparison, this study will attempt to identify only those parameter 
thresholds indigeneous to their respective sbt-minute profiles, in the hope of retaining 
as many mesoscale atmospheric features as possible. This study also proceeds with the 
presumption that by improving six-minute profiles, the hourly consensus will be im¬ 
proved. 

In order to place any credence on these thresholds it was necessary to analyze the 
largest data base possible. Unfortunately due to the combination of power outages, 
international project requirements, and equipment maintenance downtimes, the data 
base is limited. Also, in order to compare these thresholds the direction and .speed of the 
prevailing winds would have to be the same. By using the profiler consensus data, a 
10-hour period on 02 May 1990 (0600-1600 UTC) was chosen which displayed the most 
uniformity over time in the 5000-7000 m levels. Based on findings from the previous 
section, this is the highest elevation which the radar can feasibly measure reliably, i.e., 
within 3-4 m s-\ The upper atmosphere (above 7 km) was assumed to be less variable 
than the lower levels (no major synoptic-scale traveling features could be identified in the 
NMC upper air analyses), and the lower atmospheric layers below 5 km were probably 
just as variable as at any other time period which could be chosen (again, there were no 
identifiable wind shifts from frontal passages, for example, which could be identified). 
The manufacturers consensus data for this analysis period is presented in Fig. 25, and 
does in fact indicate very uniform conditions. 

The sensitivity of the receiver electronics is expected to be a limiting factor. This is 
expressed as a ratio of returned signal strength to noise. Signal-to-Noise (S/N) was first 
calculated for each of the profile gates from 

SfN = Power — [Noise -b 10 log„(a))], 

where Power is the integrated power of the peak signal, Noise is the signal strength of 
noise under the peak, and to is the bandwidth to compensate for the noise not being in¬ 
tegrated. In the analysis, (o = 256, for 256 different sampling points in the frequency 
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Fig. 25. Profiler Hourly-consensus Wind Data from 0605 to 1703 UTC 02 May 
1990: This consistent period represents the data set used to calculate 
specific thresholds from the spectral moments. 
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spectrum. Next, the U and V wind components from the six-minute files were decom¬ 
posed into the original radial velocities. Each of these radial velocities are plotted versus 
their associated S/N value in Fig. 26, all heights included. The result is a large concen¬ 
tration of points above -19 dB, and, conversely, a very broad scatter of points below -19 
dB. Statistically the standard deviation decreases dramatically, as noted in Table 3. The 
S/N thresholds in Fig. 26 are slightly different in each of the beams, e.g., -18 dB might 
be used as a value instead of -19 dB-which is presumed to be related to the direction 
of the prevailing wind, out of the north, which would favor the X-beam. Further studies 
on larger data sets of profiler wind data could better define this threshold. 


Table 3. S/N THRESHOLD STATISTICS 


Velocity 

Compo¬ 

nent 

<-21 dB 

-21to-20 dB 

-20to-19 dB 

-19to-18 dB 

>-18 dB 

X 

29.94 

23.86 

10.36 

8.84 

6.47 

Y 

41.72 

31.09 

24.45 

10.16 

10.20 


As discussed earlier, the profiler resolves horizontal winds by sampling large volumes 
of air by offset beams displaced by several meters from one another. Therefore, a 
threshold is needed to throw out data where supposedly one beam will lock on a heli¬ 
copter, airplane, or atmospheric noise and possibly even ionospheric noise. Thus, to 
relate the two beams a power difference is calculated, simply one returned power value 
subtracted from the other. Figures 27 and 28 represent the power difference between the 
X and Y-bcams, for above 5000 m and less than 3000 m respectively. Above 5000 m the 
beam power differential docs not affect the mean velocity nor the standard deviation of 
the data points, that is to say, the density of markers seem to be equally dispersed 
around the mean wind velocity component. This is presumably because the signals are 
relatively so much weaker. Below 3000 m, however, a shift in the mean velocity is evi¬ 
dent when a 5 dB differential is exceeded; the mean velocity and standard deviation are 
plotted out separately in Fig. 29. In the X-beam, the shift in mean velocity is quite ev¬ 
ident above 5 dB, but not so evident in the Y-bcam (a shift of only 1 m at 6 dB). The 
standard deviation throughout the plots, represented by the length or range of the 
shaded bars, show no apparent trend. Again, the X-beam statistics are presumed more 
prominent due to the orientation to the prevailing winds. Subjectively, by eyeballing 
individual wind barb plots which highlight suspected 'bad' winds, the 5 dB value can 
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be refined to 6 dB, i.e., more good data in relation to the consensus is retained. After 
all, an objective of this study is to make the six-minute profiles as representative and as 
useful as possible. 

Finally, the velocity variance was examined. This proved the most difficult thresh¬ 
old to determine. The premise was that the manufacturer's algorithm could not distin¬ 
guish between two different peak velocities, and, consequently, would average out the 
two peaks-efiectively masking the true wind velocity by the presence of another signal, 
such as rain or helicopters. This is illustrated by Fig. 30 at 914 to 1214 m in the spec¬ 
trum analysis. 'The vertical bar indicates the velocity which the manufacturer's algo¬ 
rithm calculated, yet there are two peaks, one of which is false and the other which yields 
the data desired. Unfortunately, there were only five sets of spectra with this feature 
which were obtained during the first month of operation, therefore only a subjective 
analysis based on a very small data set could be achieved. By comparing the velocity 
variance contained in the sLx-minute data to five sets of spectra, a velocity variance of 
3.0 m 5-' was chosen; and it was only in the lowest 2 km where the velocity variance ever 
exceeded this value, possibly due to extreme turbulence within the inversion layer. 

If each of these criteria are applied to the six-minute profiles, the six-minute profiles 
should resemble the manufacturer's consensus. Fig. 31 shows these side-by-side com¬ 
parisons of the six-minute data and the manufacturer's consensus. Each of the colors 
represent data which would be thrown out by the three thresholds determined earlier in 
this section: red barbs arc winds which exceed the velocity variance, cyan barbs exceed 
the power differential between the two beams, and green barbs exceed the minimum 
power requirements. Black wind barbs conversely are classified as the good data. It 
appears evident from this final product that quality control within the six-minute data 
is definitely useful, in that the black-colored barbs match well with the consensus data. 
Generally, the height of the radar is limited by the return strength of the signal, w'hich 
is related to the amount of humidity present. Preliminary data from Kadena airfield site 
in Japan during the Tropical Cyclone Motion experiment (TCM-90) proves that higher 
heights are obtainable in a deep, moist atmosphere, as seen in Fig. 32. There are bad 
data in the lower atmosphere as well, which are also related to beam difference problems 
and excessive turbulence. 
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Plot of Various Power Difference Thresholds Between Off-vertical 
Beams: Combined Heights from 5197-9197 m for 0605-1505 UTC 02 
May 1990. 
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riot of Beam Power Difference vs Velocity Standard Deviation and 
Mean: Combined Heights below 3000 m for 0605-1505 UTC 02 May 
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Fig. 30. E.xaniple Spectrum Depicting Velocity Shift due to Double Peak 
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Fig. 32. Representative Consensus Data from TCM-90 
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C. INVESTIGATION OF VARIABLES AFFECTING CONSENSUS 

The two variables which can be feasibly changed within the manufacturer's consen¬ 
sus procedure (i.e., software modifications) are the minimum number of points accepted 
for the consensus and the velocity threshold.2 The minimum consensus number the¬ 
oretically cannot be improved because the occurrences of incomplete hours (satellite in¬ 
terference automatic shutdown periods, power outages, frontal passages, etc.) prohibits 
reducing the value below four. This is dictated by research done by NOAA's Wave 
Propagation Laboratory (as outlined in their contract specifications to purchase a net¬ 
work of wind profilers from UNISYS). Therefore, an attempt was made to improve 
consensus by varying only the velocity threshold from 2.0 m s-\ 

Subjectively, Fig. 33 illustrates that by varying the velocity threshold from 1.0 to 
5.0 msr' there generally is little effect on the final hourly consensus profile. However, 
when the velocity threshold is diminished below 1.0 m s-\ there are noticeable deviations 
regarding direction of the individual wind barbs from the Tycho-provided consensus, 
which is most notable below 2000 m. To validate the visual conclusions, the difference 
between the new consensus (obtained wth the new threshold) and the manufacturer's 
consensus was calculated. Then the error, or "departure", of each individual six-minute 
datum was computed by subtracting the six-minute velocity from the new consensus. 
Statistically, differences between the six-minute u-velocity and v-velodity components 
and the consensus confirm what is seen graphically, in that there arc numerous "depar¬ 
tures" greater than 2.0 m s-'. Graphically, there arc also notable deviations from the 
Tycho hourly-consensus when a velocity threshold of 5.0 m s~' is exceeded, which are 
again reflected in the statistical data. For a velocity threshold above 5.0 m s~' there are 
numerous "departures" of the sbe-minute data from consensus, i.e., "departures" greater 
than 2.0 m y' 

These results indicate that there remains a possibility of reducing the velocity 
threshold below 2.0 m s"', in order to improve the reliability of the hourly data. This 
threshold was not necessarily important for the period on which this study concentrated, 
however, in other atmospheric conditions, implementation of a diminished velocity 
threshold may be useful. For example, this technique of reducing the velocity threshold 
might be more effective higher in the atmosphere, where the signal is weak and, subse¬ 
quently, the consensus number would be low. 

2 Personal conversation with Tycho Industries' technicians. 
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The data also indicates that the profiler accuracy can fall within 2.0 m !r\ since a 
consensus can be obtained with the velocity threshold set at 1.0 mr-‘. A larger data set 
of rawinsonde, or more accurate truth data, would probably confirm the accuracy of the 
profiler to be within 1-2 m 
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Fig. 33. Comparable Hourly Consensus Profiles from 02 May 1990.: DilFcrcni 
hourly consensus profiles are plotted utilizing dificrcnt velocity thresholds 
varying from 0.1 to 10.0 ms~K The hourly consensus based on the man¬ 
ufacturer's 2.0 ntsr' velocity threshold is plotted as "Radar Profiler Con¬ 
sensus". The rawinsonde data associated with this period is plotted as the 
"Sonde" profile in the left column. 
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IV. METEOROLOGICAL ANALYSIS USING SIX-MINUTE PROFILER 

DATA 


With the quality of six-minute data now established, as well as some limitations 
which constrain the accuracy, the presence of high-frequency meteorological phenomena 
can be examined. Reviewing the 30 March 1990 to 05 May 1990 period of profiler op¬ 
eration, only one distinct subperiod in the pressure, temperature, wind speed and direc¬ 
tion, and moisture data seemed to indicate the presence of any high-frequency 
fluctuations. The subperiod occurred oh 02-03 May 1990. Five-minute meteorological 
data from 03 May 1990 recorded at the Monterey Bay Aquarium Research Institute 
(MBARI) is presented in Appendix G. The wind is highly variable from 0000 to 0700 
UTC, and periodic pressure perturbations occur at 0800 and 1100 UTC in the MBARI 
data. The pressure perturbations (approximately 1 to 1 1/2 mb over one-half hour) 
though subtle, are large compared to the normally smooth diurnal change experienced 
in the Monterey Bay area. These pressure perturbations have been verified with other 
high-temporal resolution barograph data from the NFS Meteorological Laboratory. 

A. VELOCITY AND VELOCITY PERTURBATION FIELD ANALYSIS 

* 

The most commom display of six-nunute profiler data is the "time versus height' 
field of velocity. An alternative and possibly more useful way to examine the presence 
or absence of high-frequency phenomena is to examine the velocity perturbation field; 
i.c,, the vector velocities remaining after a mean value has been removed. Both of these 
methods will be utilized. 

Although the six-minute data of horizontal wind components can yield useful and 
reliable information, detection of high-frequency, low-amplitude fluctuations can be dif¬ 
ficult due to the relatively high mean horizontal wind speed. However, vertical motion 
measurements could reveal high-frequency perturbations, since the mean vertical motion 
along the West Coast under the influence of the subtropical high pressure ridge is on the 
order of centimeters per second. In fact, periodic motions reveal themselves clearly in 
Fig. 34, a time-height series of vertical motion. A dsmarcaticn line is drawn in Fig. 34 
to separate general regions of upward and downward motions. The regions are charac¬ 
terized by maximum vertical velocities on the order of 1 m s~K The demarcation line 
appears to be both oscillatory and periodic, approximately 2 hours and 10 minutes. 
Because the approximate structure of the oscillations are too small to be sampled by si- 
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multaneous measurements from the available weather network, some principles of 
time-space conversion, as presented by Forbes et al. (1989) was implemented. Forbes 
states that "the object of time-space conversion is to take frequent measurements from 
a single location while the prevailing synoptic winds [and high-frequency propagating 
phenomena] pass overhead." When time-space conversion is applied to Fig. 34, an 
internal-gravity or Kelvin-Helmholtz-like wave is suggested. It would be an error to 
assume that this high-frequency phenomena is southward propagating simply because 
the prevailing winds are out of the north. Kelvin-Helmholtz and internal-gravity waves 
traveling along a stable boundary could easily propagate against the prevailing winds. 
More information from other sites would have to be presented to determine a direction 
of wave propagation. 

Given the presence of some periodic forcing in the vertical motion, it is anticipated 
that corresponding periodic forcing is present in the horizontal velocities as well. To 
examine the periodic forcing of the horizontal velocity, velocity perturbations were cal¬ 
culated by subtracting the mean velocity from the individual six-minute observations 
(see Fig. 35). The mean was calculated at each level from the long-term consensus for 
a ten-hour period beginning 0600 UTC 02 May 1990. The horizontal velocity pertur¬ 
bation components reveal a similar periodicity as noted in the previous paragraph for 

the vertical velocity, although it is not as well defined. As expected from the "typical" 

* 

pattern of spring/carly-summer winds in the area, there was the presence of the land to 
sea breeze shift (regions highlighted by a solid line at approximately 2000 m within 
Fig. 35). However, Fig. 35 shows perturbed winds from the southwest between 
4000-6000 m from 0923 to 1012 UTC and from the northeast from 1317 to 1406 UTC, 
which correspond to two of the periods dominated by downward motions in Fig. 34. 
These regions of consistent, large perturbation velocities in the direction of the prevailing 
winds also have been highlighted in Fig. 35 by regions bordered by thick dashed lines. 
Note that the direction of perturbated winds in these latter regions arc opposite in di¬ 
rection to each other, yet both occur during corresponding times of downward motion 
in Fig. 34, which is hard to explain given the limited data and therefore left for future 
study. However, the directionality of these perturbations (aligned northeast-southwest) 
suggest that, considering conditions over the ocean are uniform and the perturbations 
are due only to the high-frequency source, the source is located to the northeast, or in 
the general direction of the pass between Atherton Peak and Fremont Peak. 
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Fig. 34. Plot of NPS Pronier Vertical Wind Component for 02 May 1990: One 
inch equals 1 m j"'. .Arrows have been olTsci from the vertical to eliminate 
overlap. 
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. 35. Plot of NPS Profiler Horizontal Wind Perturbations for 02 May 90: One 
inch equals 10 tn s'^. 















I laaa isos tsn i3<i' isos isso isst is42 issa isas 1400 i4oa i4ia i4ia i44i 1447 i4as i4aa laos isii isia iso4 isso lasa lass i64s 

^a\es from 02May90 







B. THERMAL ADVECTION FROM SIX-MINUTE DATA 

To better understand the high-frequency phenomena observed by the vertical ve¬ 
locity and horizontal perturbation velocity analysis in the previous section, some esti¬ 
mate of the high-frequency thermodynamic properties of the atmosphere would be 
useful. Although the profiler does not measure thermal properties of the atmosphere 
without a special acoustic modification (RASS), some aspects of the thermal structure 
can be inferred from profiler wind data. In particular, a single wind profiler measures 


du 

du 

dz ’ 

dt ’ 

dv 

Jl. 

dz ’ 

dt ’ 

dw 

dw 

dz ’ 

dt ' 


From these quantities, the thermal advcction can be estimated, given some rather severe 
assumptions about the (low. 

To determine the thermal advection, it is assumed that all Coriolis, earth curvature, 
frictional, and vertical acceleration terms arc negligible compared with those involving 
the vertical pressure force and the force of gravity. This leads to the hydrostatic 
equation presented by Holton, 1979 


where clz is the layer thickness, dp the pressure difference between layers, p the density, 
and g the acceleration of gravity. The thermodynamic state of the atmosphere at any 
point is determined by the values of pressure (p), virtual temperature {Ty), and density 
(p) at that point. If the fluid can be regarded as incompressible these field variables can 
be related to each other by the equation of state for an ideal gas as presented by 
Campbell (1973) 

P — P^d^Vt ( 2 ) 

where Rj is the gas constant for dry air. Solving Eq. (2) for p , and then substituting 
Eq. (2) into Eq. (1) 


54 





dz 


RdTy dp 
g P ' 


( 3 ) 


Integrating Eq. (3) throughout the layer gives 


fv. =f 


dp 

g P 


( 4 ) 


The term on the left side of Eq. (4) is a definite integral, and Rj and g are constants; 
therefore, from Eq. (4) 





( 5 ) 


Next, assume that temperature is linear throughout the layer and Eq. (5) becomes 



( 6 ) 


Two more equations are needed, Holton presents the hypsometric equation, which re¬ 
lates pressure differences to height or thickness changes as 


Az = 


(Alnp)/?^ _ 
g Iv 


( 7 ) 


Finally, an assumption is made that the accuracy of the geostrophic wind as an ap¬ 
proximation to the observed wind is a good one, and the difference is negligible, i.e., 
Coriolis force exactly balances the pressure gradient force. Gill (1982) presents the 
geostrophic approximation in component form as 

= (8) 


g dz 
f dx' 


( 9 ) 


or in the vector form, 
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(10) 


The vector difTerence is clearly 

aP, = (11) 

which by substitution of Eq. (10) into Eq. (11) yields 

^ r X V;^2a - 2,) = J A X V;^A2). (12) 

Substituting Eq. (7) into Eq. (12) 

= (13) 

but, RJgl ln(y7,//>j)] = A 2 / 7 ', by Eq. (7), so Eq. (13) now becomes 

= xv„r^. (M) 

By definition of the thermal wind by Holton (1972), and substitution of Eq, (14) into the 
definition 



Finally, temperature advection can be defined as 


iL 

di 


rAd. = -y-^T, or 
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( 18 ) 


^ dT dT 

^Adv “ 3 ^ V . 


Substitute Eq. (16) and Eq. (17) into Eq. (18), and replace u with (the consensus ve¬ 
locity) yields 



f^v 
Az g 


]-v,C- 


fTy 
Az S 


]. 


(19) 


Realistically, and are not known so u and v, the profiler-derived winds, are 
substituted. Ty was calculated from a formula derived from Stull (1988), 
ry = (1-f-.bl^)^, using data from the 1300 UTC 03 May 1990 rawinsonde launch. 
Large deviations c way from geostrophic balance are possible for six-minute time scales 
so the geostrophic assumption is likely to produce substantial errors in V»VT. Note 
that this analysis will highlight regions of strong vertical shear in the horizontal 
ageostrophic (perturbation) winds. 

When temperature advection values are calculated from Eq.(19) with the sbe-minute 
profiler data, the result is generally nothing short of a random distribution of minimum 
and maximum values. How’ever, throughout a two-hour period on 03 May 1990, de¬ 
picted in Fig. 36, a startling consistent pattern of concentrated minimum and maximum 
values lying along sloped lines is revealed. Again, the profiler does not measure tem¬ 
perature, but if it is assumed for a moment that Fig. 36 truly is temperature advection, 
then it would be expected to also have packed isotherms in the regions of the large 
advection values (i.e., packed isotherms would be needed in order that homogeneous 
atmospheric vertical motion would carry "drastic" changes in temperature.) The high 
values of temperature advection in Fig. 36, both positive and negative, are probably 
associated with tighter temperature gradients, but not directly. The slight pressure per¬ 
turbations that arc induced by wave activity also set-up temperature gradients, similar 
to the packed isotherms Gossard (1990) revealed regarding turbulence in a stratified gas 
(sec Fig. 37), These isotherm gradients in turn set up w'eak velocity perturbations which 
are subscqucnily detccicd by ihe profiler. Although this phenomenon is easily depicted 
using this temperature advection computation, it is actually only an artifact created by 
the associated weak velocity perturbations observed by the profiler. 

Subsequently, the vertical velocity for the period represented in Fig. 36 was also 
plotted, as Fig. 38. Lines indicating regions of the highest density of extreme temper¬ 
ature advection values from Fig. 36 were transferred to Fig. 38. There is a tendency for 
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Temperature Advection Field: Reveals areas of high vertical shear In the 
horizontal winds. 











Fig. 37. Packing of Isotherms due to Turbulence in a Stratified Gas: (Source: 

Gossard (1990).) 

vertical convergent winds (Fig. 38) to correspond to the large values of temperature 
advcction (Fig. 36), which supports the idea that substantial changes in the horizontal 
velocity perturbations arc occurring across these layers. 

C. SYNOPTIC OVERVIEW FOR 02-03 MAY 1990 

Synoptic conditions at the 500 mb and surface, representative of the period 02-03 
May 1990, are presented in Fig. 39 and 40, respectively. The combination of an upper- 
level low height center over southern Arizona and a high pressure ridge along the West 
coast force prevailing 500 mb winds out of the north over California. Surface conditions 
along the West coast arc influenced by a surface pressure trough extending from the 
Arizona-New Mexico border to northern California. This trough creates north by 
north-westerly flow principally parallel to the California coast; however, southerlies ly¬ 
ing along the California central coast to the south of Monterey are against the synoptic 
scale pressure pattern. Thi.s suggests a possible mesoscalc eddy or low pressure center 
to the north and west of Monterey. 

Local hourly-consensus data from the profiler for the period 02-03 May 1990 are 
presented in Fig. 41 through Fig. 44. Winds out of the north prevail from 2000-6000 
m. The land-sea breeze shifts are evident in the lowest 2000 m of the atmosphere. 




















Fig. 38. Vertical Velocity Field: The heavy lines represent lines of large temper¬ 
ature advection value concentrations transposed from Fig. 36. 




Specifically, westerlies depict the sea breeze crossing the Marina shoreline of the 
Monterey Bay, which persist from approximately 2100 UTC to 0300 UTC. From 0400 
UTC to 2000 UTC, the land breeze sets up from the southeast, very much aligned to 
topography of the Salinas River valley. Hourly observations at Salinas, California (20 
km down into the valley, southeast of the Monterey Bay) are presented in Table 4. The 
onset of the land breeze is evident at 0500-0700 UTC on both days in the hourly surface 
observations, denoted by the large speed decrease in the northerlies. 

The sbc-minute data for this period of land/sea breeze shift is presented in Fig. 45. 
As seen from the sbc-minute data presented in Fig. 45, there is a "burst" of easterlies in 
the lowest 4000 m after about 0414 UTC 02 May90 (as well as on 03 May 90 at ap¬ 
proximately the same time). (This "burst" is more evident in the hourly consensus data. 
Fig. 41.) The easterlies are really a vector summation of two different circulation pat¬ 
terns: (1) the land breeze from the south by southeast, and (2) the prevailing 
northeasterly synoptic-scale surface flow. This "burst" of easterlies occurs because the 
sea breeze (out of the west) dams the synoptic northeasierlies between Atherton and 
Fremont Peaks, and when the sea breeze collapses the northcasterlies now flow through 
the pass easily and abruptly. (Refer to Fig. 2 for topography of the local area.) 

D. A PROBABLE EXPLANATION OF OBSERVED HIGH-FREQUENCY 
PHENOMENA 

Although the determination of the source and propagation nature of the oscillatory 
wave seen in Fig. 34 arc beyond the scope of this study, a hypothesis is required to plan 
future research work. Two possible explanations are that some type of internal-gravity 
or Kclvin-Helmholtz instability waves exist. Although the wave pattern certainly re¬ 
sembles the propagation nature of an internal-gravity wave, a calculation of the 
Vaisala-Brunt frequency at 2000 to 4000 m yields a value or9.7 x 10-’j-'-which equates 
to an oscillatory period of approximately two minutes. This period differs considerably 
from the two hour period which is seen Fig. 34. Therefore, 1 conclude that the high- 
frequency phenomena is probably due to shear instability, specifically from the shear 
induced by the reversal of the sea breeze (i.e., land breeze southerlies under the prevail¬ 
ing northerlies aloft). This is particularly convincing as the oscillations begin shortly 
after the land breeze begins. Additional six-minute profiler data can be used to confirm 
this hypothesis in the future. 


61 

















Fig. 41. NFS Radar ProFiIer Hourly-consensus Data for 02 May 1990 
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Fig. 42 . NFS Radar Profiler Hourly-consensus Data for 02 May 1990 
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Fig. 43. NFS Radar Profiler Hourly-consensus Data fur 03 May 1990 
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Fig. 44. NFS Radar Profiler Hourly-consensus Data for 03 May 1990 
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Table 4. SALINAS HOURLY SURFACE WEATHER OBSERVATIONS 
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V. CONCLUSIONS 


This study determined the accuracy of the NPS profiler to be 3-4 ms-\ i.e., the 
standard deviation values taken from the atmospheric layers which exceeded a .90 
corelation coefficient in Chapter III. However, by obtaining a consensus with a 1 msr^ 
velocity threshold, which was similar in all aspects to the manufacturer's consensus, 
suggests that the accuracy of the profiler is significantly better than 3-4 m llius there 
is an indication that the profiler does meet its contractual requirements. This disparity 
is easily explained by (I) the dislocation of sensors, (2) the recent evidence ('Thomson 
and Williams 1990) that rawinsondes are less accurate than radar profilers, and (3) the 
extremely small data set available for this study. 

Reliable data was found within the six-minute profiles, it is not as random as it first 
appears. The parameter thresholds determined in Chapter IV are very effective for 
identifying bad data and can be useful in improving the individual sbt-minute profiles. 
The thresholds also can be used individually to isolate certain physical characteristics 
of the air column such as turbulence. And when all the data is assimilated, a good pic¬ 
ture of what is occurring in the PBL (up to 6 km) is obtained. Ironically, even the 
presence of bad data has a meaning, such as the presence of a dry, subsiding air mass. 
It is then not enough to say the data is bad, it must be understood why the data is bad. 

Although the data quality checks which were developed to clarify the six-minute 
data were very successful, the presumption that they wou^^ also improve the hourly 
consensus was not proven with as much success; the hourly consensus was not dramat¬ 
ically affected by the improvement of six-minute data. 

Despite the large vertical spatial resolution (250 to 1000 m) in the radar data, the 
overall radar performance appears to be very good. A great deal of detail can be drawn 
from the wind profiles. It has been showTi that high frequency meteorological wave 
phenomena can be detected with the six-minute data, but other fluctuations in the local 
meteorology can also be detected. The six-minute data can accurately depict the dra¬ 
matic and sudden changes in the land/sea breeze circulation. Furthermore, there are 
complex topographic interactions which may be understood by analyzing the six-minute 
data. The high-temporal resolution and accuracy of the profiler in determining three 
dimensional meteorology is certainly an exciting addition to the NPS Meteorology cur¬ 
riculum. 
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VI. 


RECOMMENDATIONS 


This study has been able to demonstrate the quality and usefulness of the profiler 
data during its early period of operation. However, there is considerable potential for 
future and more complete meteorological studies. The following list of recommen¬ 
dations are things that may help check the quality of the data and better define the de¬ 
tailed meteorology sampled by the profiler. 

A. MEASUREMENT VALIDATION AND DATA QUALITY IMPROVEMENT 

• Improve the manufacturer's algorithm for spectral averaging based on local mete¬ 
orology. 

• Allot a rawinsonde system dedicated to work with the profiler. 

• As Thomson (1990) indicates, the uncertainty in profiler-derived \vTnd speeds and 
direction is clearly significantly less than the natural variability of the •"'ind. Rather 
than using rawinsondcs which are highly prone to sampling errors, another radar 
profiler, LIDAR, or SODAR should be used to accurately gage NPS profiler 
performance-not the rawinsonde system. 

• Investigate techniques such as those listed in Appendix D which can improve data 
quality outside the hourly-profiles. 

B. ACADEMIC RESEARCH 

• Use a 50 MHz profiler for comparison with the NPS profiler. It should detect 
winds higher in the atmosphere. 

• Investigate the possibility of acquiring a complete network of at least three 
profilers, either through outright purchasing or coordination with other agencies 
or academic institutions. For example, San Jose State University might be inter- 
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ested in purchasing a profiler, placing it within the Monterey Bay local, and then 
recover the data remotely via telephone lines. 

• Add the NFS 2 kHz SODAR to the operational data base. 

• Coordinate 5 kHz SODAR operations with the Physics Dept. 

• Pursue implementation of the manufacturer's software regarding Q, calculations. 

• Investigate the possibility of determining a Power to Humidity algorithm. 

• Rigorously employ the Radio-Acoustic Sounding System (RASS) modification 
package within the profiler. 

• Make all remote connections to automate delivery and storage of the profiler data, 
especially the capability to upload data into the NPS IDEA Laboratory. 

• Get more students involved. Draw real-time information into the classroom. 

• Conduct a detailed study of land/sca breeze and topography interaction along the 
Monterey Bay. 

• Incorporate real-time ground- and space-based remote sensing data, such as wind 
profiles from the NPS Doppler wind profiler, in as many of the meteorology 
courses as is possible. 
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APPENDIX A. EQUIPMENT CONFIGURATION 
A. WIND PROFILER 



Fig. A - I. Functional Diagram of Wind Profiler 


'I'hc wind profiler consists of a two-stage transmitter, antenna, receiver, 
transmit/rcccivc (T/R) switch, and processor with its firmware. Fig A-1 illustrates the 
basic operation of a Tycho-engineered wind profiler. On command from the Real-Time 
Processor (RTP) the antenna controller sets up the proper phasing for the desired an¬ 
tenna beam position (i.e., vertical, 15® southeast, or 15® southwest). The T/R pulse sets 
up the 'f/R switch for transmit mode. The RTP sends a prompt signal (or pulse) to the 
rcceivcr/modulator, which produces a pulsed RF signal. This is then disconnected from 
the antenna during the transmit cycle to prevent overloading. The system is then set to 
the "receive" configuration. The receiver is reconnected and the T/R switch directs the 
weak return signal from the antenna to the receiver. The receiver amplifies the signal 
and extracts the in-phase and quadrature (or sine and cosine) phase components from 
which the radial velocity is derived. Filtered outputs are sent to the RTP, which per¬ 
forms the various processing steps necessary to produce the radial wind components, 
noise levels, signal power strength, and velocity variance. The Profiler Data Handler 
(PDH) is a VAX 3100 smart terminal which processes these parameters to derive time- 
height wind profiles, Q profiles, wind speed and direction profiles, radial velocity vari¬ 
ance versus height, and signal-to-noisc (S/N) profiles. The Keyboard Display Unit 


73 














(KDU) has four functions: (I) it provides the ability to edit operating parameters, (2) 
displays individual mode spectra, (3) monitors profiler operation, and (4) controls 
profiler operation. 

B. RAWINSONDE 

The upper air profiling system utilized for this analysis integrated the 
Vaisala-DigiCORA receiver-processor with Vaisala RS80-series rawinsondes. The spe¬ 
cific software version installed within the DigiCORA system was the Corex version. 
Measurements are based on the use of a frce-ascent balloon rawinsonde. The 
rawinsonde, or "sonde”, transmits data to the station at 402 MHz. Refer to Fig. A-2. 
Pressure, temperature, and humidity are measured by sensors in the sonde. Each of the 
capacitive sensors controls the frequency of an oscillator through an electronic switch. 
Relationship between the frequency and corresponding parameter value is established in 
the calibration process prior to launch. (Specifications of the individual sensors are 
provided in Fig. A-3. The location of the sonde is determined by means of a third gen¬ 
eration Navaid (Omega + Very Low Frequency (VLF))-based windfinding program 
called "multifrequcncy composite VLF". In the composite solution more than one fre¬ 
quency from each Omega transmitter and one frequency of the USSR VLF network are 
used. Navaid signals are relayed to the ground station for processing and for wind pro¬ 
file computation. Briefly, the signals arc processed in a correlation processor to derive 
their relative phases from which the sonde's movements can be calculated. Wind speed 
and direction over the range of 0 to 180 ms-' is computed with an accuracy of 0.5 
ms-', regardless of balloon distance from the ground station. 






Fig. A-2. Vaisala DIgiCORA Upper Air Sounding System: (Source: Vaisala 


DigiCORA Repair Manual) 
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Table A-1. RAWINSONDE METEOROLOGICAL SENSOR TECHNI 


CAL DATA (from Vaisala Corporation) 


Pnssure 


Type 

Capacitive Aneroid 

Measuring Range 

1060 hPa to 3 hPa (mb) 

Resolution 

0.1 hPa 

Accurao' 

0.5 hPa 

Temperature 

Sensor Type 

Capacitive Bead 

Measuring Range 

+ 60‘C to -90*C 

Resolution 

o.rc 

Accuracy 

0.2*C 

Humidity 


Sensor Type 

HUMICAP* Thin Film Capacitor 

Measuring Range 

OtojOOVoR// 

Resolution 

\%RH 

Lag 

1 s (6 m J-‘ flow at 1000 hPa, +20“C) 

Accuracy 

2% RH 








APPENDIX B. DERIVATION OF THE RADAR EQUATION FOR CLEAR 

AIR PROFILERS 


The radar equation relates the received power (echo power) to radar parameters and 
target parameters. If the radiating source radiates isotropically, that is, equally in all 
directions, the power density (P,, power per unit volume) at range R from the radiator 
is 


Pr = 



(B-l) 


where P, is the peak transmitted power. Radar antennas do not radiate isotropically 
however, but rather transmit the energy in a beam. The energy density actually received 
at any given range and direction relative to isotropic, is called the gain (G) or 

G = , (5_2) 

Msotropic 

If the target scatters isotropically, the power intercepted by the receiver is 

Pr “ (iA -X ^ 


where o is the back-scatter cross section of the target, and A, is the effective antenna 
aperature, and a is the transmitter antenna cfilcicncy. 

The widths of the main beam in the north-south and east-west directions, 6 and (f> 
respectively, are generally specified by the w-idths at the half-power points. These widths 
depend mostly upon the size of the antenna and wavelength of the radiation, although 
the shape of the antenna and power distribution across the antenna also influence these 
widths. To a good approximation 

sin d = --, {B - 4a) 

sin 4> = * 7 ^. {B - ^b) 






Fig. B- 1. Illustratiun of Gconietrj' for Beam Volume Calculations 


where L^ and Lj arc the linear dimensions of the antenna in the directions that corre¬ 
spond to 0 and (f ), and X is the transnutted wavelength. For narrow beams, the sine of 
the angle can be replaced by the angle itself: 

0 == , and . {B - 5a) 


= (B-Sb) 

•’2 

■fhe width of the beam between the first nulls (the nulls between the main beam and 
sidclobcs on either side of the main beam) is about twice these half-widths. The solid 
angle corresponding to these half-widths (Q) using liqs. (B-5a) and (B-5b) is 


The area of a rectangular planar array (/t^) by substitution into Hq. (B-h) is the same 
as the physical area: 


/f = L,/-2 = 



(/f- 7) 
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But if the power distribution across the antenna is tapered, the effective area will be 
slightly different from this. In general, the effective aperature is smaller than the phys¬ 
ical aperature, the ratio being called the antenna efficiency (a): 

a = ,or Ag = ctA. {B — 8) 

The antenna gain can be related to the aperature by making the simplifying approxi¬ 
mation that the actual beam paiiem can be replaced by a rectangular beam pattern that 
has a given gain G over a width equal to the half-width and zero everywhere else. In this 
case 


Pt 



4nR^ 


but by substituting Eq. (B-7) into (B-9) 



Substituting Eq. (B-10) into Eq. (B-3) returned power becomes 


{B-9) 


(5-10) 


p, = p,G,A^fl, ~ f . (B- H) 

IvTT I\f I\f 

but, a, = a, = a, 5, = /?, = R, C7, = (7, = G, and A, = A, = A., because the antenna is the 
same hardware used for the transmitter and receiver, therefore 


TbTF- 

For a distributed target such as encountered in meteorology, we must sum over all the 
scatterers in the sample volume. This sample volume is defined by the radar range re¬ 
solution, AR = ct/2 (where t is the pulse length, and c is the propagation speed of radio 
energy), and by the half-widths of the radar beam 0 and </>: 

2 2 
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Assume the scatters don't "shadow" each other and that they don't have appreciable 
absorption. Then, the radar reflectivity is the sum of all the as in the volume divided 
by the volume: 


E” 

)7 = -y-; orcr = ?/K (5-14) 

By substitution of Eq. (B-14) into Eq. (B-12) the radar equation becomes 

Utilizing Eq. (B-10), by substitution into Eq. (B-15) 

Utilizing Eq. (B-13), Eq. (B-16) becomes 

crR^X^ 

P, = P,A^ (x^rj = P,Aa^cr —. (5-17) 

16;r^5^ 87r5^ 

The above assumption can be made more accurate by using a more appropriate form 
of the gain equation. If a Gaussian beam shape is used, instead of a rectangular shape, 
the fact that the edges of the sample volume arc not illuminated as powerfully as the 
central portion can be compensated by adding a shape factor: 

5, = iP,A VcT —~ ] X iShapeFactorl (5-18) 

For a planar array according to Gossard and Strauch (1983) the shape factor is 
4I9{X^IA). Substituting the shape factor into Eq. (B-18) the radar equation evolves into 

P, = P,Aa^cx —(5-19) 
187r5^ 

Eq. (B-19) is appropriate for the scattered signal. The radar is, of course, also sensitive 
to noise. Noise comes from a variety of sources, the main ones being system noise (also 
called Johnson noise) and cosmic noise. System noise can be controlled somewhat by 
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the design of the radar but it can never be completely eliminated. Noise power is typi¬ 
cally expressed in terms of the equivalent temperature that a blackbody would emit at 
the given frequency to account for the power observed. The noise power (?;,) is 

P^=kBiaTc+Ts), (5-20) 

where k is Boltzmann's constant, B is the receiver bandwidth (which is designed to match 
the transmitted bandwidth, approximately 1/t at the half-power points), Tc is the cosmic 
noise temperature, and is the system noise (also called the receiver noise). The factor 
a arises because inefTiciencics in the system affect the noise power received as well as the 
signal power: 


P,= P,A,a'^cx^ 


n _ 

UR^k[aTc+Ts) ‘ 


(5-21) 


This is a basic equation used by the Tycho profiler and is consistent with research done 
by Rottger and Larsen (1990). Finally, the radar constants are taken out of the equation 
and grouped together for ease, as used in the body of the study: 

“ %nk{aTc + Ts) . (5-22) 


where is a "radar constant". 
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APPENDIX C. RADAR WIND PROFILER TRIGONOMETRY 


ZENITH 



Fig. C- 1. liorizont.'il Velocity Calculation from tlie Three Radial Velocities 


Figure C-l illustrates how the horizontal wind components can be derived from the 
profiler-measured radial velocities. Fig. C-l illustrates three i. 'ams which are utilized in 
either the low- or high-power mode. The "zenith", or vertical beam, measures the w- 
component of the wind directly. The radial velocities of the other two off-vertical beams, 
Xj^„ and have not only a vertical component, but also horizontal components 
related by the following equations: 

y — 3' ' 4- y " V — V ’ ^ V “ 

~~ ^RAD ^ ^FjiO ‘ RAD ~ ‘ RAD ^ * RAD 

Xrad == — (- F;, sin 02 -1- w cos Oj) = — (w cos 0, + u„ sin 0,) 
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Note: The negative sign compensates for the scheme depicting velocities toward the 
radar as positive. If w, and (the three radial velocities processed by the 
profiler) are known; then 

«/? = (- YRad - vi” cos 0,) / shn 0,; and 
v/j = - (- - w cos O 2 ) I sin 02 

Note: A negative sign must be applied to v^j, to compensate for the usual horizontal 
cartesian coordinate system, i.e., v to the north. 0, and 0^ are the angles which the off- 
vertical beams are offset from the vertical. These angles are currently set at 15°. «« is 
the Eastward wind component computed from the radial velocities, and is 
Northward wind component computed form the profiler radial velocities. 
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I/*/- upXcos(11.5*) + v,,x sin(II.5‘) 

Component 1 Component 2 

Component 3 Component 4 

Fig. C-2. Trigonometry Relating Coordinate Axes in the Horizontal Plane 
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Fig. C- 3. Trigonometry of Axis Translation for Site Correction 


Site Alignment 
Error i tl.5® 
Site Alignment 
I 55" 


Nots: Correction for site alignment to files before 03 May 1990 2200 (Ul'C) must be 
applied. After 03 May 1990 2200 (UTC), the elevation and orientation arc corrcct--no 
correction needs to be applied. 
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APPENDIX D. CONSENSUS TECHNIQUES 

Because we know that spurious data can enter the profiler system, much more 
readily that a direct-measuring system, automated quality control (QC) procedures have 
been developed. This appendbc contains excerpts from Brewster (1989). These proce¬ 
dures take advantage of the consistency of the winds in time and height (on the scale 
of sample separation, 250 meters in height and 6 minutes or 1 hour in time). The QC 
algorithms for the profiler data calculate an expected value to the measured value. The 
expected value is typicaliy a function of the measurements surrounding the datum in 
time and space. It may also depend on some sort of "first guess" provided by a forecast 
model or expressed as a function of the surrounding data. If the difference between the 
expected and measured values exceeds a threshold, then the datum is considered suspect. 
Various procedures may be followed to determine if the datum arises from a plausible 
but extreme meteorological event or if the large difference is caused by a bad datum 
among those used to calculate the expected value. In this appendix several techniques 
that can be applied to wind profiler data arc described-first that which is developed for 
use by the Tycho profiler, then some others that can be applied in-house by forecasters 
and are included for reference as techniques to apply in future studies using the NFS 
profiler. 

For wind profiler data distributed to data networks in the future the QC must be 
done in real time, that is, the data for the current hour must be checked before data from 
the next hour arrive. Estimating measured values is more difficult because the data are 
unevenly distributed in time-plentiful data from past times, but none from future times. 
For some applications, more precise "two-sided" quality control, done after future data 
have arrived, is not only possible but desirable. Two factors merit consideration: 1) how 
long after the data collection can one wait before running the QC, and 2) how sensitive 
is the application to data errors. 

A. CONSENSUS AVERAGING 

The 1-hour wind averages produced at NFS are derived from the 6-minute radial 
velocity samples by a technique known as consensus averaging. Frimarily an averaging 
process, this technique also has built-in QC. As suggested by the name, not all the data 
are averaged, but only those that mutually agree within a threshold. A 2.0 m s** window 
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is used for the velocity. We review the effect of this technique by considering these ten 
6-minute radial velocity samples: 

3.3 4.6 20.9 2.8 3.6 4.1 -6.8 3.4 4.0 22.2 (ms“’) 

The average of all samples is 6.2 m s"'. Looking at the indi/idual points you can see that 
6.2 m s-' does not at all represent the typical sample value during the hour; it has been 
severely contaminated by three outliers. Now if we collect with each sample all other 
samples agreeing within ± 1.5 m s*', that is falling within a 3.0 m s-' window, we get the 
groups from Table D-I. The greatest number of samples in a group is seven, a tie among 
five groups. The group of seven coming latest in the hour (indicated by the arrow) is 
then selected for the consensus; in this case all five groups of seven are identical so the 
position within the hour doesn't matter. The consensus average is the average of that 
group, 3.7 m s~\ which almost anyone would agree is representative of the entire sample 
set. 

Table D-I. FUNCTIONAL CONSENSUS PROCESS: Example of individ- 
_ ual sample cross-comparisons involved with consensus process. 


sample 

samples within window 

3.3: 

2.8 

3.3 

3.4 

3.6 

4.0 

4.1 

4.6 

4.6: 

3.3 

3.4 

3.6 

4.0 

4.1 

4.6 


20.9: 

20.9 

22.2 






2.8: 

2.8 

3.3 

3.4 

3.6 

4.0 

4.1 


3.6: 

2.8 

3.3 

3.4 

3.6 

4.0 

4.1 

4.6 

4.1: 

2.8 

3.3 

3.4 

3.G 

4.0 

4.1 

4.6 

-6.8: 

-6.8 







3.4: 

2.8 

3.3 

3.4 

3.6 

4.0 

4,1 

4.6 

4.0: 

2.8 

3.3 

3.4 

3.6 

4.0 

4 1 

4.6 -<— 

22.2: 

20.9 

22.2 







In practice, mathematical shortcuts speed calculation of the average value that has 
the most samples lying within the threshold window. If no group has at least four 
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members to form a consensus (as might be found if the data are made up of random 
'-’alues), then the average value is reported as zero. 

Consensus averaging handles isolated spikes in the data quite well, as in the exam¬ 
ple. It also handles cases in which the returned atmospheric signal is too weak to .be 
detected. Then, the data samples are nearly random, and a missing value will be re¬ 
ported. 

B. COMBINED MEDIAN FILTER AND SHEAR CHECK 

As the name suggests, this routine consists of two checks applied sequentially to 
hourly averaged data. In the median check (see Fig. D-1) data are gathered from adja¬ 
cent hours (the current hour and the last two available hours) and adjacent levels. The 
median (the reported value of which half the samples are greater and half are less) of this 
collection of data is then computed for the two horizontal wind components. If the 
difference between the datum and median is greater than a threshold, the median is re¬ 
calculated using only those data from the current hour and the last hour. This is done 
to preserve data in the event of a strong wind shift with time as might be found in a 
dramatic trough passage. Data are flagged bad if the observed datum and the recalcu¬ 
lated median also differ by more than the threshold. The threshold value depends on the 
height of the datum (greater thresholds at greater altitudes), the wind speed (greater 
thresholds for high speeds), and the time difference between the datum and ihe data 
collected for the median computation (greater threshold when there are missing data in 
previous hours). In the current implementation, the median check is quite loose; it is 
designed to flag the most widely varying data but may pass along a substantial number 
of spurious data. The philosophy is to discard the implausible values before the shear 
check and to leave the more difficult decisions to that second step. 

The vertical consistency check uses a gatc-to-gate shear threshold, which is a func¬ 
tion of height, distance between the observations, wind speed, and the diflcrence in wind 
direction (greater .shear is allowed if the wind direction is unchanging). A typical 
threshold is about 20 ms-' per kilometer or about 13 kts per thousand feet. 

If the magnitude of the vector shear between the first and second gates exceeds the 
threshold, the shear between the second and third gates is checked. If the second shear 
is nearly the same magnitude as the first and in the opposite direction, the data in the 
second gate are flagged bad (both horizontal components). This situation is illustrated 
in Fig. D-2. The magnitude of the shear between the first and second gates, 22.6 kt, is 
excessive, thus prompting a check on the shear between the second and third gates (27.5 
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Fig. D — 1. Data Points from Mhich Median is calculated in the Median Check 

ki). The two shears are comparable in magnitude and nearly opposite in direction (the 
components have opposite signs). Thus the light northwest wind in the second gate is 
flagged. 

If the wind in the third gate tends to continue the shear between the first and second 
gates, a straight line is plotted through the u- and v-components of winds in the first and 
fourth gates. If both components in the second and third gates lie reasonably close to 
their respective lines, they arc accepted as valid measurements. This is the condition 
shown in Fig. D-2b, where the sample winds imply strong shear across a frontal 
boundary. On the other hand, if cither component in the second and third gates lies far 
from the line, both winds are flagged as bad. Fig. D-2c illustrates the case of two spu¬ 
rious winds between two valid winds. The line for the u-components lies along the 
horizontal axis. Although the middle two u-components lie close to the line for u, the 
middle two v-components lie far from the line for v. Thus the middle two winds are 
flagged. Note that the scale for u and v in Figs. D-2b and D-2c has been held constant 
to facilitate comparisons. 

'fhe practical result of these tests is that a layer of winds containing four or five 
measurements that pass the median check and corroborate each other will be retained 
even if there is substantial shear at the layer interfaces. Moreover, two adjacent, spuri¬ 
ous winds supporting each other are likely to be flagged correctly as bad, as in the above 
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Fig. D-2. Tiiree examples of tlie vertical consistency check.; (a) A bad wind 
bciween two good winds, (b) good winds with strong shear, and (c) bad 
winds in gates 2 and 3. 
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example; three such winds may not be. Conversely, if a layer of anomalous but valid 
winds is only samples by two or three gates, those wind observations may be erroneously 
rejected. 

The vertical consistency check is applied once to the low mode data and then to the 
entire profile, high and low modes meshed together. Checking in each case begins with 
the lowest gate and proceeds upward. Special processing with linear extrapolation of 
neighboring data is used at the lowest and highest gates in the profile, should large shear 
be indicated there. 

To re-emphasize, the most powerful check in the system is the vertical consistency 
check. Substantial changes from one hour to the next are generally allowed to pass 
through the median check but are rejected later if there is insufficient continuity in the 
w'ind shift with height. 

C. OPTIMAL INTERPOLATION 

Statistical objective analysis, also known as Optimal Interpolation (01), is widely 
used for analysis of meteorological variables on quasi-horizontal surfaces. In fact, 01 
is used today at many national weather centers, including the National Meteorological 
Center (NMC) and Fleet Numerical Oceanography Center (FNOC), for model initial 
analysis. This method utilizes statistical properties of the observations and the back¬ 
ground field (also called a "first guess"). The background is usually a model forecast but 
climatological averages or persistence can serve the same purpose. The 01 analysis takes 
into account the positions of the observations relative to each other and relative to the 
analysis point; it has an extrapolative property that aids analyses near data edges. 

Optimal interpolation can be used for quality control. The datum under scrutiny is 
compared with a value analyzed at the same location by means of 01. Surrounding 
observations (but not the central one) and the background field contribute to the ana¬ 
lyzed value. If the difference between the analyzed and observed value exceeds some 
threshold, the datum requires further testing. Since the 01 method provides an estimate 
of the analysis error, the threshold can be specified as a function of the analysis error. 
This allows for more variation when observations arc sparse or distant from the datum 
in question and less variation when the observations arc dense or located very close to 
the datum. 

Large discrepancies between the observed and analj'zed value may arise either be¬ 
cause the subject obser%'ation is erroneous or because one of the observations contrib¬ 
uting to the analysis is erroneous. To find out which is true, the data are re-analyzed 
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several times, removing one of the surrounding observations each time to see if that 
produces agreement (Fig. 3). The removed observation is returned for subsequent re¬ 
analyses so that only one observation is absent for each re-analysis. If no agreement is 
found for any re-analysis, the datum being examined is flagged. Should the removal-of 
one of the observations produce agreement, the removed observation is flagged as sus¬ 
pect and is not used in any subsequent QC analysis. Suspect flags are internal to the 
QC program and are generally not part of the QC output. Typically, a datum is flagged 
only when it is the subject datum. 

This technique just described can be used to check data in a time-height cross- 
section as well as on the horizontal surfaces for which it was originally designed. The 
correlation statistics must then be specified in time-height space. 01 time-height quality 
control of profiler data in real time presents a special problem. Because there Is an at¬ 
tempt to always examine the latest hour, the analyses are not centered in time (i.e., it is 
always analyzing on the edge of the time-height domain). A second, often vexing 
problem is that spurious winds are not always isolated-repeated errors in a given beam, 
extending over several gates, can cause the rc-analysis technique to fail. 

The majority of cases can be handled properly if data selection for each analysis is 
done to minimize the possibility of multiple bad datj, points being included. An 01 QC 
algorithm developed at PROFS makes two passes through the data with the flagging 
threshold decreasing on the second pass. In the first pass, data are checked in sequence 
from the lowest gate to the highest. The two nearest data points arc collected from each 
of the numbered sectors shown in Fig. D-4, except those which have not yet been 
checked (those in the shaded area-Scctor I and pan of the sector containing the subject 
point P). The eight data points which correlate best with the point being exanained (as 
determined from historical data) are used in the analysis. The flagging system is the 
same as described in the paragraph above. The second pass proceeds as the first, but 
data are accepted from all the sectors. In this pass it is presumed that if there are any 
spurious data above the point being examined (i.c., not yet checked on this pass), they 
have been flagged bad in the first pass or flagged suspect when they w'erc used as a 
contributing observation on the second pass. 

D. RECURSIVE FILTER 

Quality control can be built into the analysis itself by adjusting the observation 
weights based on the observation diflTcrence from a first quess or, more commonly, from 
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Fig. D- 


Orioinal analysis for observation A 


CD 


CD 

CD ^ 
CD CD 


A= observation being checked 
analysis location 

1...5 neighboring 
observations 


First reanaivsis 



CD CD 


Analysis is redone at point A 
using observations 2-5. 
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Fig. D-'l. Data Search Strategy for Optimal Interpolation Analysis (01) 
QC: Sector searching assures an even distribution of data around the da¬ 
tum being checked. Data within shaded area arc excluded in the first pass; 
all sectors arc searched in the second pass. 

the prc\ious analysis in a multiple pass (successive correction) anal)sis method. In this 
approach, data points arc not discarded but arc given a very small weight. 

The recursive filter analysis uses bilinear interpolation to place the obscr\cd data 
(actuallv the difrcrcncc from a first guess field) onto a rcgularl) spaced grid. To smooth 
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the analysis, a numerical filter (similar to a weighted running average) is applied to the 
grid-point values. Then the process is repeated using the last analysis as a first guess. 
With each pass, the filter is modified so that less smoothing is done. Quality control is 
introduced by weighting each observation according to its difierence from the previous 
analysis. One can discard an observation if its "quality weight" becomes very small, or, 
if the data density is high enough, one can presume that observations with very small 
weight wDl be overwhelmed by those with large weight (good quality observations). 
Because the recursive filter requires far less computation than optimal interpolation and 
no data searching it is much faster than 01. The speed depends upon the number of 
re-analysis passes necessary and the number of filter passes performed for each analysis 
pass. A good first guess can make the number of analysis passes very small (e.g., five) 
and can greatly increase the accuracy of the quality control decisions. 

This technique has been tested on six-minute profiler data at PROFS. The analysis 
is done on a lOO-m by 6-minutc grid and the first guess consists of the median of 11 
observations around each grid point. A time-to-space conversion which equates six 
minutes to 100 meters is used, and data are discarded (QC flag set to "bad") if their 
weight is less that lO'h 
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APPENDIX E. INDIVIDUAL PROFILER COMPARISONS AND 
RAWINSONDE LAUNCHINGS 
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Fig. E- I. Profiler-RaHinsonde Wind Coinparisoiis for 20 April 1990 
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Fig. E-2. R.'iu insonde Soimding from 20 April 1990 
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E“3. Profiler-RaMinsonde Wind Comparisons for 16U() UTC 23 April 1990 
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Launch time ; 90 04 23 
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E-4. Ra>vinson(le Sounding from 1600 UTC 23 April 1990 
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Fig. E-7. Profiler-Raw insonde Wind Comparisons for 2(K)0 UTC 02 May 1990 
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E-8. Ra>un.soiide Sounding from 1900 UTC 02 M-iy 1990 
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Fig. E-9. Prolller-Rauinsomle Wind Comparisons for 0100 UTC 03 May 1990 
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E- 10. RaMinsonde Sounding from 0100 UTC 03 May 1990 
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li>ijiiiiiiiuiuiuliiiiiiiiiuiiiuiiiIiiiiiuiiuiiimiiluiuiuuiuiuiiLiuuuuiiiiiniiliuiiuiiuuiiuiiliuiiiiuiiiiuiiuliuuuuuiuujiJiiuuuiiiuuiiu] 


O 

O 

O 

CO 


O 

O 

o 


o 

o 

o 


o 

o 

o 


o 

o 

o 

o 


o 

o 

o 

oo 


o 

o 

o 

lO 


o 

o 

o 


o 

o 

o 

rvi 


(sj9|suj) epnimv 


Pig. E - II, Profiler-RaMinsoiuIe Wind Comparisons for 0700 UTC 03 May 1990 
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Fig. E- 12. Ra^unsonde Sounding from 0700 UTC 03 May 1990 
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Fig. E- 13. Profilcr-Rauiiisomle Wind Coinp.nrisons for 13(10 UTC 03 M.iy 1990 
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E- [‘j. RaHiiisonde .Sounding from 13UI) UTC 03 May 1990 
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Profiler — Rawinsonde Comparison Winds 

Launch iime; 90 05 05 at 18 26 GMT. Site Lat: 36.250, Long: 123.23 
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Fig. E- 15. Profiler-Ra«insonde Wind Coinp.msuns for 19U0 UTC 03 May 1990 
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Fig. E- 16. Rawla.sonde Sounding Trom 1900 UTC 03 May 1990 
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APPENDIX F, NFS PROFILER SITE SURVEY RECORD AND LAYOUT 


Naval Pootgi'aduate School 
Depat'ttiieril; of Ocearionraphy 

I 

* 19 OPfl 90 




Frofiii CoiiUitancJot* Kui't Schncbsln, NOAD (Coda UC/Sn> 

Toi Profuuuor UilliaM Shaw (Code MIO 

Sub)i nZIHUlM t POSITION OF FORT ORD PROFIt.ER 

End I (1) nzitiuith t P'oaition Dotartiiinat ion of Foi't Drd Profllai* 
(ii> Fort Ord Profiler Survey Uorltaheuta 
(3) Pouitions t Oeacriptions of Control Stationa 


I. At the r{?(|ueBt of LCOR Dan Strecd, USN; Nr. Jariina Cherry and 
1 coMpletud a aurvey to datarMtne a peodutic aziiiiulh ar>d poult iori 
foi* the profiler array at Fort Ord. The ruaulta are aumiiiarisod 
below and dlocuuned in Encloauru (i). 


OslMtith of Side 336* 07.0* (linear error SSX a 3.3* > 
(wuuturn aide) • 

Asihiutti of North 021* 07.0’ (auautiiing array in aquare) 
South Dlaponal 


Array Center 
Poult ion on 
NAD 03 datuhi 


N 36® 41’ A3.716" (drculai« error 95v!) 
W 121® AS* 29.oaA" (a 3 cm relative to) 
(otation DAIRY. I 


Elevation SI Miitera above ooa level (NIIVD 29) 

Array Top 


2. Note that the elevation waa deterMined froM only one point 
(iitatlon DAIRY). CachinQ an independent cliecR froM another 
ourvey point, it lu beat to auuuMo that elevation iu only n<^'C'd to 
±1 Muter. Doth the pouition and azlMutti were deterMined uBtn,'| 
redundant MeauvireMunta in a Ieaut-aquareu adJuatMunt, which pave 
tlio prucieion eutiMatuu cited above. 

3. Note that we aot a now aurvey MOnuMont (DAIRY OFFSET) on the 
til I 1 ovo.*Iool<ltia the profiler. If needed, the axiMiitli could tie 
chochod froM there in the future. 






Kurt J. Schnebele 
CoMMandor, NOAA 
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MARIN 8 to AIRPORT BEACON 

rRlTZCIIE AAP PT URD 


Fig. F-1. Atmospheric Profiler, Fort Ord Survey Observ.ntiou Sketch: ('on- 
ductcd 19 April 1990 by Cdr. Schncbclc (no scale). 
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(deg C) (degrees) 


yVPPENDIX G. HIGH-FREQUENCY METEOROLOGICAL DATA 



Fig. C-I. Recorded Data for 03 May 1990.: From the Monterey Hay Aquarium 
Research Institute. 
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